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SUMMARY 


This  document  describes  the  construction  of  methods 
for  predicting  the  hydrodynamic  characteristics  of  sub- 
mersible vehicles  in  pitching  attitudes  with  undeflected 
tailfins.  Following  an  approach  used  successfully  in  mis- 
sile aerodynamics,  a set  of  models  was  built  and  tested  to 
obtain  systematic  data  over  relevant  ranges  of  geometry  and 
flow  conditions.  These  data  were  then  used,  in  conjunction 
with  theoretical  results,  to  develop  the  prediction  methods. 
The  methods  deal  with  the  characteristics  of  individual 
vehicle  components  (bodies,  tails)  and  with  their  mutual 
interactions  when  combined  into  complete  configurations. 

Each  method  is  presented  in  self-contained  form,  along  with 
directions  for  its  use.  It  is  shown  that  the  methods  pro- 
vide good  accuracy  for  prediction  of  hydrodynamic  character- 
istics . 
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B1 

B2 

B3 

Cy 


NOMENCLATURE 

coefficient  of  zero-order  terra  in  expression 
for  nondimensional  load  or  moment 

coefficient  of  first-order  term  in  expression 
for  nondimensional  load  or  moment 

coefficient  of  second-order  term  in  expression 
for  nondimensional  load  or  moment 

fin  aspect  ratio  - two  panels 

total  span  of  two  fin  panels  placed  together 

distance  between  vortices  or  span  of  vortex 
wake 

conical  base  section,  1 caliber  in  length 
conical  base  section,  2 iCalibers  in  length 

conical  base  section,  3 calibers  in  length 

» 

fin  local  chord  length 
crossflow  drag  coefficient 

turbulent  flat  plate  skin-friction  coefficient 
based  on  exposed  surface  area,  Sw 

rolling  moment  coefficient 

base  drag  coefficient 

skin-friction  drag  coefficient,  based  on  SR 
axial-force  coefficient 
lift  coefficient 
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NOMENCLATURE  - CONTINUED 


'M 


'N 


T 

CA 


CB 


CC 

Cl 

C2 

C3 

C4 

d 


f 

fi (Cq/sJ 

f2  (A) 
f3  <lLE/d) 


pitching-moment  coefficient 
normal-force  coefficient 
fin  root  chord 


fin  tip  chord 

cylindrical  center  section,  2-1/2  calibers  in 
length 

cylindrical  center  section,  1-1/2  calibers  in 
length 

cylindrical  center  section,  3 calibers  in 
length 

cylindrical  center  section,  2-1/2  calibers  in 
length,  consisting  of  CA 

cylindrical  center  section,  4 calibers  in 
length,  consisting  of  CA  + CB 

cylindrical  center  section,  5-1/2  calibers  in 
length,  consisting  of  CA  + CC 


cylindrical  center  section,  7 calibers  in 
length,  consisting  of  CA  + C3  + CC 


maximum  diameter  of  body  cross  section;  also 
reference  length 

diameter  of  effective  or  equivalent  base  (i.e. 
diameter  of  base  section  of  body  where  the 
axial  flow  separates) 


transverse  force  per  unit  length 


multiplicative 
effect  of  fin 

factor  used 
span 

in 

B(T) 

for 

multiplicative 
effect  of  fin 

factor  used 
taper  ratio 

in 

^N 

B(T) 

for 

multiplicative 
effect  of  fin 

factor  used 
position 

in 

CN 

NB(T) 

for 
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^(B) 


NOMENCLATURE  - CONTINUED 

empirical  factor  used  in  C.. 

MB 

ratio  of  normal  force  curve  slopes:  tail  on 
body/tail  alone 

apparent  mass  factor  of  Munk 

length 

length  of  body  exposed  to  potential  flow  (i.e., 
distance  from  nose  to  xaa) 

axial  distance  from  start  of  base  section  to 
fin  root  chord  leading  edge 

normal  force 

hemispherical  nose 

ellipsoidal  nose,  1 caliber  in  length 
ellipsoidal  nose,  2 calibers  in  length 
dynamic  pressure,  P^vf/2 
local  body  radius 
maximum  radius  of  body 
body  radius  at  fin  leading  edge 
body  radius  at  fin  trailing  edge 
ratio  of  normal  forces:  tail  on  body/tail  alone 
CnT(B)/CnT 

constant  term  of  tail  normal-force  amplifica- 
tion factor 

first  order  term  of  tail  normal-force  amplifica- 
tion factor 
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NOMENCLATURE  - CONTINUED 


Re  Reynolds  number  based  on  body  length 

ReCR  Reynolds  number  based  on  root  chord  of  fin 

Red  Reynolds  number  based  on  body  maximum  cross- 

sectional  diameter 

Re^  crossflow  Reynolds  number,  Vjdsina/v 

Rega  Reynolds  number  based  on  length  of  body 

exposed  to  potential  flow 

RM  rolling  moment 

s fin  maximum  semispan  measured  from  body 

centerline 


s 


local  cross-sectional  area  of  body 


s 


b 


w 


Vr 


X 

X 


cross-sectional  area  of  equivalent  base, 
irdJ/4 

area  of  body  planform 

planform  area  of  portion  of  body  in  separated 
crossflow  (i.e.,  from  x_„  to  end  of  body) 

BC 

2 

maximum  area  of  body  croRB  section,  rrd  /4 , 
reference  area 

fin  planform  area,  single  panel 

surface  area  of  body 

body  volume 

free-stream  velocity 

axial  or  chordwise  coordinate 

axial  distance  to  center  of  pressure 

axial  location  of  moment  reference  center 

,1 

axial  location,  measured  from  beginning  of  body  j 

base  section,  where  the  axial  flow  separates  | 

I 
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NOMENCLATURE  - CONTINUED 


y 

Y 

a 

au 

r 

& 

a 

Ax 

n 

a 

p» 


axial  distance  on  body,  measured  from  nose, 
where  crossflow  begins  to  separate 

center  of  pressure  of  body  carryover  loading, 
^N 

NB(T) 

lateral  co-ordinate 

lateral  distance  to  center  of  pressure 
angle  of  attack 

angle  of  attack  of  upwash  flow 
vortex  strength 
tail  deflection 

incremental  contribution;  used  in 
development  ' 

correction  to  8^  to  give  ^T(B) 

correction  for  three-dimensional  effects  on 
body  crossflow  drag 

fin  taper  ratio 

free-stream  density 

complex  distance  between  a vortex  and  its 
image  inside  a cylinder 


Subscripts 

A aft  tail  location 

At  incremental  contribution  due  to  fin  aspect 

ratio  effects;  used  in  development 

b base 

B body 
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BT 
B (T) 


c 


C 

d/b 

P 

1 

LE 

M 

nl 

N 

pot 

Re 

T 

T (B) 
o 

a 

X 


NOMENCLATURE  - CONCLUDED 


complete  body-tail  combination 

carryover  on  the  body,  in  the  presence  of  the 
tails 

centroid  of  body  region  between  fins 
center  section 

contribution  attributed  to  span  effects?  used 
in  development 

forebody,  or  forward  tail  location 

incremental  contribution  due  to  location  of  fin 
on  base  section;  used  in  RT(Q)  development 

denotes  value  at  leading  edge  of  root  chord 

mid  tail  location 

non-linear  term 


nose 

potential  term 

incremental  contribution  due  to  Reynolds  number 
effects?  used  in  development 

tail  alone 

tail,  in  the  presence  of  the  body 

denotes  value  at  zero  angle  of  attack 

denotes  differentiation  with  respect  to  a 

incremental  contribution  due  to  fin  taper  ratio 
effects?  used  in  development 


(Reverse  Page  8 Blank) 
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I.  INTRODUCTION 


A recurring  problem  in  submersible  vehicle  design  has 
been  the  lack  of  accurate  methods  for  predicting  configura- 
tion hydrodynamic  characteristics. 

Without  accurate  methods,  submersible  design  is  an  un- 
certain process,  requiring  extensive  testing  and  modifica- 
tion of  models  and  full-scale  vehicles  before  an  acceptable 
configuration  is  defined.  As  a result,  development  times 
are  lengthened  and  both  time  and  cost  schedules  may  be 
adversely  affected.  It  cannot  be  said  that  the  availabil- 
ity of  accurate  methods  is  a sufficient  condition  for 
removal  of  all  obstacles  to  effective  designs.  However, 
it  is  certainly  a necessary  one. 

The  main  obstacles  to  construction  of  such  methods  have 
been  (a)  the  inability  of  classical  hydrodynamic  theory  to 
handle  viscous-flow  situations  and  (b)  the  lack  of  suitable 
data  to  guide  semi-empirical-method  construction. 

Classical  theory  has  been  proven  quite  reliable  in 
cases  where  configuration  flow  fields  are  not  significantly 
affected  by  viscous-flow  phenomena,  i.e.,  by  boundary  layers 
and  wakes.  For  submersible  vehicles,  whose  base  regions 
produce  thick  or  separated  boundary  layers  and  their  associ- 
ated wakes,  and  which  deploy  their  control  fins  in  precisely 
those  viscosity-dominated  regions,  classical  theory  provides 
only  a rough  guide  to  performance.  Of  course,  theoretical 
results  may  be  supplemented  and  modified  through  reference 
to  experimental  data.  However,  if  the  data  are  unsystematic 
(as  are  those  from  specific  vehicle  designs)  there  is  little 
chance  of  constructing  methods  having  wide  utility  for  broad 
classes  of  vehicles.  What  is  required  is  a theoretical  base, 
supplemented  by  systematic  data  and  it  is  this  approach  which 
has  been  used  in  the  present  work. 

Following  an  approach  used  successfully  in  missile  aero- 
dynamic design,  Ref.  1,  the  work  described  here  Bet  out  first 
to  design,  build  and  test  a large  number  of  models  whose 
geometric  parameters  systematically  covered  the  ranges  typ- 
ical of  submersible  designs.  Data  from  these  tests  were 
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then  used  to  supplement  theoretical  results  in  the  con- 
struction of  semi-empirical  methods  for  predicting  submers- 
ible hydrodynamic  characteristics. 

To  facilitate  data  generation  and  instrumentation  de- 
sign, the  tests  were  carried  out  in  a wind  tunnel,  not  in 
water.  Following  the  wind  tunnel  tests,  however,  a limited 
investigation  was  made  in  a small  water  tunnel  to  determine 
whether  test  medium  had  an  effect  on  the  data.  The  results 
of  these  tests  are  given  in  a separate  report. 

The  contents  of  this  report  are  as  follows!  Descrip- 
tions are  given  of  the  modele,  tests,  data  generated  and 
method  construction.  Each  individual  method  is  presented  in 
self-contained  form  and  includes:  discussion  of  method  con- 
struction, description  of  use,  a numerical  example,  charts, 
graphs  and  tables,  and  a comparison  between  method  predic- 
tions and  experimental  data.  Discussions  of  configuration 
hydrodynamic  behavior  are  given  to  supplement  the  methods 
descriptions. 
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II.  MODELS 


The  models  were  designed  to  cover,  systematically, 
those  ranges  of  geometry  typical  of  submersible  vehicles, 
viz. , 

• Bodies  having  overall  slenderness  ratios  varying 

from  four  to  twelve. 

4 Noses  having  mainly  ellipsoidal  shapes,  of  varying 
length  from  0.5  to  2.0  diameters. 

• Bases  of  conical  form  (idealized  shapes)  having 

slenderness  ratios  from  one  to  three. 

' Tails  of  varying  aspect  ratio  (0.5  to  2.0),  taper 
ratio  (0  to  1.0)  and  body  diameter/exposed  tailspan 
ratio  (1.0  to  1.8) . 

More  detailed  descriptions  of  the  models  follows i 

Model  centerbody  diameter  was  7 inches.  Nose  shapes 
are  as  shown  in  figure  1.  Three  ellipsoidal  noseB  are 
available  having  length/diameter  ratios  0.5,  1.0  and  2.0, 
plus  a torpedo-like  nose  having  a profile  of  aquation  shown. 
The  noses  are  designed  to  be  readily  interchangeable.  Each 
can  be  fixed  to  the  model  centerbody. 

Model  centerbody  is  composed  of  three  separate  sections 
as  shown  in  figure  2.  The  section  CA  accommodates  the 
main  recording  balance  during  test.  Additional  sections  CB 
and  CC  can  be  deployed  forward  and  aft  of  CA  respectively. 
These  center  sections  permit  the  cylindrical  portion  slender- 
ness ratio  to  be  varied  from  2.5  to  7.0. 

A typical  model  base  is  shown  in  figure  3.  In  all, 
three  bases  were  built,  lengths  1,  2 and  3 calibers.  Each 
is  of  conical  shape.  The  bases  are  made  of  aluminum  and 
are  machined  to  accommodate  the  three-component  balances  to 
which  the  tails  are  fixed.  When  no  tails  are  required,  the 
holes  are  filled  by  dummy  plugs. 
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FIGURE  2.  CENTER  SECTIONS. 
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Tail  geometry  and  si2es  are  shown  in  figure  4.  The 
tails  are  deployable  on  the  bases  and  also  on  a reflection 
plane  which  .is  shown  in  figure  5.  This  plane  permits  test- 
ing of  half-tails  and  yields  results  which  are  representa- 
tive of  complete,  two-panel  tails.  The  tails  are  mounted 
on  three-component  strain  gauge  balances  which  were  special- 
ly-built for  this  program.  These  balances  are  shown  in 
figure  6.  They  measure  tail  normal  force,  pitching  moment 
and  rolling  moment.  On  the  reflection  plane,  tail  angle  of 
attack  is  changed  by  rotating  the  plane.  In  the  model,  tail 
deflection  is  changed  by  means  of  * ndax  plates  which  mate  to 
the  balances.  Each  tail  may  be  deflected  to  the  angles  0, 
±10,  ±20,  ±30  degrees.  During  these  initial  tests,  the 
tails  were  undeflected. 

Photographs  of  various  model  pieces  are  shown  in 
figures  7 through  10. 


FIGURE  4.  TAIL  GEOMETRY. 


FIGURE  7.  MODEL  NOSES 
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FIGURE  10.  MODEL  P 
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III.  WIND  TUNNEL  TESTS 


The  teste  were  conducted  in  the  12-Foot  Pressure  Wind 
Tunnel  at  NASA/Ames  Research  Center.  Moffett  Field,  Cali- 
fornia. Wind  tunnel  testing  was  chosen  for  the  following 

reasons . 

• Theoretically  there  is  no  difference  between  dynamic 

characteristics  in  water  and  in  air  so  long  as  the 

test  Reynolds  numbers  are  equal. 

• Testing  in  air  is  far  simpler  than  in  water,  partic- 

ularly where  electrical  apparatus  is  involved. 

A summary  of  the  test  configurations  is  given  in  Table 
1.  In  general,  each  configuration  was  tested  at  the  free 
stream  Reynolds  numbers,  4x10*  and  B.5xl0  per  foot.  Body 
and  Body  Tail  configurations  were  mounted  on  a strut;  tails 
were  mounted  on  a reflection  plane.  Tail  angles  were  varied 
from  -3®  to  +30®  on  the  reflection  plane.  On  the  body,  the 
tails  were  undsflected.  Body  angle  was  varied  from  -3®  to 
+15®.  Typical  teat  configurations  including  the  reflection 
plane,  an  isolated  body  and  a body  tail  combination  are  shown 
in  figures  11  through  13.  Before  proceeding  to  method  con- 
struction, thj»  data  were  corrected  for  zero  shifts  where 
necessary. 

The  data  generated  included : 

• Six-components  of  forces  and  moments  from  a Task 

balance  located  inside  the  body. 

• Three-components  of  force  and  moment  from  each  tail. 

In  addition,  oil-flow  visualizations  were  made  of 
several  configurations.  Still  photographs  and  videotape 
recordings  were  made. 
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T5M 

T14 


M 


M 


T2 

T9M 

T13 

Til, 

T2A 

T12, 

*10, 

T8A 

Til, 


M 


TABLE  I. 

CONFIGURATIONS  TESTED 

Til 

\ 

T14 

T8 

T9 

V 

T13 

/ 

T15 

T12 

T10 

T7 

/ 

N2 

Cl 

B2 

\ 

N3 

C2 

B2 

N1 

C2 

B2 

N2 

C2 

B2 

k 

N2 

C2 

B3 

N2 

C4 

B2 

N3 

C4 

B3 

N1 

Cl 

Bl 

N2 

C2 

Bl 

/ 

N4 

C2 

B3 

N2 

1 

C2 

1 

B2 

1 

TllM  1 

Isolated 


Isolated 


v Body  tall 
• combinations 


*< 

■ * 


NCSC  TM- 238-78 


TABLE  I.  CONCLUDED. 


B2 

1 

T12 

T 

B1 

T12 

T15 

B2 

T12 

B3 

Tllj 

T2M 

If 

t5m 

B2 

T12 

B2 

Til, 

4 


Body  tail 
combinations 


A - Aft  tail  location. 

M - Mid  tail  location. 

F - Forward  tail  location. 
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IV.  FORMULATIONS 


A 


k 


Following  the  practice  widespread  in  aerospaoe  work, 
the  hydrodynamic  characteristics  of  complete  configurations 
were  formulated  by  considering  first  the  characteristics  of 
isolated  components  (bodies,  tails)  and  then  combining 
these,  accounting  for  their  mutual  interactions. 


From  the  reflection  plane  tests,  the  following  quanti- 
ties  are  obtained. 


C„  Isolated  tail  normal-force  coefficient 

nt 

XT  Isolated  tail  ohordwise  center  of  pressure 
Isolated  body  tests  yielded 


C..  Body  normal-force  coefficient 

NB 

CM  Body  pitching-moment  coefficient 

B 

C.  Body  axial-force  coefficient 

ab 

For  body/tail  configurations,  more  complex  quantities 
are  obtained,  as  described  in  the  following. 


The  various  quantities  associated  with  the  hydrodynamic 
pitching  behavior  of  a cruciform  submersible  are  written 
below.  The  horizontal  tail  half-panels  are  treated  as 
though  they  had  arbitrary  deflection  6 at  body  pitch  angle 
a.  (During  the  tests  6 was  kept  equal  to  zero).  The 
deflection  6 is  assumed  positive  leading-edge-up  relative 
to  the  tail  neutral  position.  The  body  is  assumed  fixed  in 
roll  altitude  and  the  vertical  tails  are  undeflected,  thus 
generating  no  force. 


C 


XBT 


Ji  V™*!001 4‘ 


?T(S)i  8t 


(1) 
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Here,  R^bj  and  ^t(B)  are  the  normal“force  amplifica- 

i i . u 


tion  and  center  of  pressure  location  on  the  i 
tail  panel. 


th 


horizontal 


'N, 


RT(B)i 


T(B) 


(2) 


'N„ 


This  is  analogous  to  the  KT/B)  factor  of  reference  2. 
Hence,  RT,0)  and  ^T(B1  ar®  functions  of  a,  6.,  and 
geometry,  'i  ' *i  “ 1 

With  these  quantities  available,  it  iB  possible  to 
write  the  axial-force  coefficient  bb 


BT 


il;i[VT<B!i 


sin  ff  + CAi 


B 


(3) 


Note  that  the  effects  of  tail  leading-edge  suction  have 
been  ignored.  For  these  tails  they  are  negligible. 

The  othar  configuration  coefficients  are: 

2 - i S„ 


ch  - cm  + £ fcM  R™  yB\  oos  6,  + CM 

NBT  NB  i-.lL  NT  T(B)Ji  1 SR  N 


‘M, 


BT 


Cmb  + 


liiV1 


(B) 


T(B} 


cos  5. 


B(T) 

T_ 

R 


(4) 


+ C. 


NB(T) 


(5) 


where  C 


M is  the  carryover  body  normal  force  (related 

nB(T) 


to  of  ref.  2)  and  3?^  is  its  center  of  pressure. 

The  value  %p(B)  is  the  chordwise  center  of  pressure  loca- 
tion of  the  i^  horizontal  tail. 

To  summarize  then,  the  methods  to  be  constructed  deal 
with  the  following  quantities. 
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C.  body  axial-force  coefficient 

aB 

CN  body  normal-force  coefficient 

B 

C„  tail-alone  normal-force  coefficient 

N_ 

T 

CK  carryover  loading  on  body  due  to  presence  of 

"B (T)  the  horizontal  tails 


(B)  i 
*N 

^P(B)i 

YT(B)i 

C«B 


amplification  of  tail-alone  normal  force  in 
presence  of  the  body 

aenter  of  pressure  of  CM 

nB(T) 

chordwise  center  of  pressure  of  tail  on  body 
spanwise  center  of  pressure  of  tail  on  body 
body  pitching-moment  coefficient 
chordwise  center  of  pressure  of  tail 


Methods  were  constructed  for  all  these  quantities. 

Each  method  contains  a background  discussion , a description 
of  method  construction,  directions  for  use,  numerical 
example,  comparison  with  independent  data,  references, 
charts,  graphs  and  tables. 


(Reverse  Page  38  Blank) 


37 


V*  i 


NCSC  TM-238-78 


V.  ISOLATED  BODIES 


V.l  GENERAL  CONSIDERATIONS 

When  a slender,  cylindrical  body  is  placed  at  an  angle 
of  attack  a in  an  air  stream  of  velocity  VM,  the  flow 
may  be  considered  as  consisting  of  two  components.  One  com- 
ponent is  parallel  to  the  body  axis  and  is  termed  axial 
flow  Vaocos  a.  The  other  component  is  normal  to  the  body 
axis  and  is  the  crossflow,  V^sin  a.  At  angles  of  attack 
greater  than  about  6°,  the  boundary  layers  associated  with 
the  crossflow  separate  on  either  side  of  the  body  and  form 
a lee-side  wake  as  depicted  in  figure  14  taken  from  reference 
3.  This  wake  takes  the  form  of  a pair  of  symmetrically  dis- 
posed, counter-rotating  vortices  and  is  fed  by  vorticity  shed 
from  the  separating  boundary  layer.  Initially,  the  separated 
region  is  quite  small  and  is  located  near  the  rear  of  the 
body.  Eventually,  however,  the  separation  region  extends 
over  virtually  the  entire  body  as  angle  increases.  At  angles 
of  attack  greater  than  about  25®  an  asymmetric  vortex  flow 
pattern  may  appear.  This  report  does  not  consider  the  lat- 
ter type  of  flow  pattern. 

The  loadB  on  a cylindrical  body  at  angle  of  attack  are 
usually  described  in  terms  of  a normal  force,  an  axial  force, 
and  a pitching  moment.  In  general,  all  three  quantities  are 
linear  with  angle  of  attack  for  small  angles.  At  larger 
angles,  however,  these  quantities  demonstrate  nonlinear  be- 
havior, duoi  to  the  separated  crossflow  discussed  in  the  pre- 
vious paragraph.  The  axial  force  shows  relatively  little 
variation  with  angle  of  attack  up  to  at  least  15®  but  the 
normal  force  and  pitching  moment  can  show  quite  large  varia- 
tions, depending  upon  the  particular  configuration.  Past 
studies  (e.g.,  ref.  4)  have  been  quite  successful  in  applying 
potential  flow  theory  to  predict  the  linear  behavior  of  these 
quantities  and  a crossflow  drag  analogy  to  predict  the  non- 
linear behavior.  According  to  thiB  latter  concept,  each 
section  along  the  body  experiences  a cross  force  equal  to  the 
drag  force  the  section  would  experience  with  the  axis  of 
revolution  normal  to  a stream  moving  at  a velocity  V.sin  ot 
(ref.  4)  . 
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The  use  of  potential  flow  theory  in  the  prediction  of 
the  transverse  force  on  a body  of  revolution  is  well  estab- 
lished, having  been  initiated  by  Munk  in  connection  with 
his  work  on  airships  (ref.  5) . Munk  showed  that  the  cross 
force  per  unit  length  can  be  expressed  as 

f - q ||  sin  2a  (6) 

Thus,  applied  to  bodies  that  are  closed  at  the  base,  an  in 
the  present  test,  potential  flow  theory  yields  a zero  net 
force,  but  a nonzero  pitching  moment.  This  means  that  the 
load  on  such  a body  is  in  the  form  of  a pure  couple.  The 
center  of  pressure  of  such  a load  is  at  upstream  infinity. 
This  result  is  true  only  for  angles  of  attack  near  zero.  In 
reality,  the  boundary  layer  along  the  body  separates  some- 
where upstream  of  the  end  of  the  body.  The  location  of  this 
separation  point,  x , can  thus  be  used  to  determine  the 
amount  of  potential8Iift  on  the  body.  The  essential  linear- 
ity of  the  transverse  force  with  a at  small  angles  shown 
by  equation  (6)  is  borne  out  by  experiment.  Hence,  at  small 
angles  of  attack,  the  normal  force  on  an  isolated  body  is 
due  to  potential  lift  which  extends  back  to  the  location 


At  higher  angles  of  attack  the  crossflow-drag  analogy 
is  used  to  explain  the  nonlinear  viscous  crossflow  effects 
on  the  loads.  The  crossflow  normal  force  has  the  form 


C„  - nCd  sinso  (7) 

nl  c R 

It  is  important  to  note  two  points  associated  with  the 
viscous  crossflow.  First,  the  crossflow  does  not  begin  to 
separate  at  the  nose  but  at  some  downstream  position  x 
which  is  a function  of  the  angle  of  attack.  Second,  vortex 
size  and  strength,  in  general,  increase  toward  the  rear  of 
the  body.  Thus,  the  proper  distribution  of  crossflow  force, 
expressed  in  the  term  Cd  (x)  is  necessary  to  properly 

c 

account  for  the  nonlinear  pitching  moment. 

In  the  following  sections  these  concepts  will  be  uti- 
lized to  aid  in  developing  methods  for  the  calculation  of 
body-alone  normal-force  coefficient,  axial-force  coeffi- 
cient, and  pitching-moment  coefficient. 
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The  results  of  the  flow  visualization  investigation  of 
the  isolated  bodies  brought  out  an  interesting  feature  of 
the  flow  over  the  one  caliber  base.  Figures  15  and  16  are 
oil  flow  photos  of  the  one  caliber  base  and  reveal  the 
presence  of  two  concentrated  symmetric  vortices  located 
beneath  the  base.  The  presence  of  this  vortex  structure 
can  entirely  dominate  the  hydrodynamic  characteristics  of 
bodies  with  this  size  base.  Even  at  low  angles  of  attack 
this  highly  viscous  region  means  that  linear  potential  flow 
theories  do  not  provide  accurate  load  estimates.  The  loca- 
tions for  the  horizontal  fins  are  seen  clearly  in  Figure  15 
and  the  direction  of  the  flow  at  that  location  means  the 
fins  will  be  in  a region  of  separated  flow,  greatly  reducing 
their  effectiveness.  This  is  discussed  in  more  detail  in 
Section  VII. 1.  When  applying  the  various  methods  to  bodies 
with  a one  caliber  base,  the  presence  of  this  wake  structure 
should  be  borne  in  mind.  A similar  investigation  of  longer 
bases  did  not  reveal  the  presence  of  any  concentrated  vortex 
structure  - the  flow  was  smooth  and  attached  out  to  the  end 
of  the  base. 


V. 2 BODY-ALONE  NORMAL-FORCE  COEFFICIENT 

This  method  permits  calculation  of  the  normal-force 
coefficient  on  isolated  bodies  in  incompressible  flow  for 
angles  of  attack  up  to  15°.  The  method  consists  of  two 
steps.  The  first  determines  the  normal-force-curve  slope 
from  the  effective  base  and  potential  flow  theory.  The 
second  calculates  the  nonlinear  viscous  force.  Comparisons 
between  the  method  and  experimental  data  show  good  agree- 
ment. 


The  normal  force  on  an  isolated  slender  body  of  revolu- 
tion at  angle  of  attack  is  made  up  of  contributions  due  to 
potential  lift  and  viscous  crossflow  effects.  At  small 
angles  of  attack  the  normal  force  is  well-approximated  by 
potential  flow  theory  and  is  essentially  linear  with  angle 
of  attack.  When  the  body  exceeds  an  angle  of  attack  of  about 
t>  a lee-side  vortex  wake  forms  and  the  normal  force  becomes 
nonlinear.  Past  work  haB  shown  (e.g.,  ref.  4)  that  the 
total  normal  force  can  be  quite  successfully  described  by  a 
combination  of  potential  flow  theory  (to  describe  the  linear 
behavior)  and  a second  nonlinear  term  derived  from  the  cross- 
flow-drag  analogy  described  in  the  previous  section.  Hence, 
the  normal-force  coefficient  is  assumed  to  be  of  the  form 

CN  ■ axa  + a2a*  (8) 
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FIGURE  15.  OIL  FLOW  PHOTOGRAPH  OF  BASE  Bl 
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FIGURE  16.  OIL  FLOW  PHOTOGRAPH  OF  BASE  Bl. 
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The  following  boundary  condition  has  been  applied: 

CM  (a  - 0)  - 0 (9) 

NB 

Therefore,  taking  the  derivative  of  equation  (8) , the  result 
is  derived  that  is  equal  to  the  zero  angle  normal-force- 

curve  slope,  i.e., 


The  coefficient  a2  is  the  value  given  by  determining  the 
crossflow  drag  of  the  particular  body. 

Description  of  Method 

Potential  Force.  Por  a body  not  closed  at  the  rear, 
the  normal  force  can  be  expressed  as  (Ref.  6) 

NP  ■ q £ d£  sin  2a  “ 2qas£  (11) 

where  Sb  and  db  are  the  cross-sectional  area  and  diameter 
respectively  of  the  base.  In  this  relation  the  normal-force 
coefficient  becomes 


C 


N 


pot 


(12) 


where  SR  is  the  maximum  cross-sectional  area  of  the  body. 
Stated  in  terms  of  the  normal-force-curve  slope,  potential 
flow  theory  thus  gives 


in  the  configurations  used  in  the  present  test  the  body  is 
closed  at  the  bass.  This  implies  that  d^  and  hence  Cm 
are  zero.  However,  the  experimental  values  of  CNfl  are 

non-zero  and  equation  (13)  can  be  used  to  determine  the 
diameter  of  the  "effective  base"  caused  by  boundary  layer 
separation  or  thickening.  Then,  from  the  body  surfaoe 
coordinates,  the  location  of  axial-flow  separation  oan  be 


45 


NCSC  TM-238-78 


i 

I 


i 


i 


determined.  This  result  is  useful  when  tail  fins  are  placed 
in  this  region  as  it  will  be  possible  to  estimate  how  muoh 
of  the  fin  is  immersed  in  separated  flow  and  how  much  is 
immersed  in  the  external  potential  flow. 


Examination  of  the  data  from  the  present  test  revealed 
that  the  normal- force-curve  slope  was  a function  of  the  base 
length  (or  equivalently  dS/dx  of  the  base)  and  of  the  fore- 
body slenderness  ratio  (that  portion  of  the  body  ahead  of 
the  conical  base  section) . Correlation  curves  were  prepared 
from  the  experimental  values  of  , along  with  equation  (13) , 


to  determine  the  location  of  x#a,  the  position  of  boundary- 
layer  separation  of  the  axial  flew.  These  curves  are  given 
in  figure  17.  The  distance  x_a  is  measured  from  the  start 
of  the  conical  base  section.  Hopkins,  in  a previous  study 
(ref.  7)  was  able  to  correlate  xBa  measured  from  the  nose 
with  the  distance  from  the  nose  to  the  point  at  which  dS/dx 
has  a maximum  negative  value  (xt).  This  is  essentially 
similar  to  the  curves  in  figure  17.  All  the  bodies  used  in 
Hopkins'  study  had  smooth  and  continuous  variation*  in  the 
cross-sectional  area.  In  the  bodies  examined  in  the  present 
experimental  study,  the  location  of  the  maximum  negative 
value  of  dS/dx  always  occurs  at  the  start  of  the  base 
section.  Thus,  lp/d,  the  forebody  slenderness  ratio,  corre- 
sponds to  x,  in  Hopkins'  work.  However,  sinae  models  of 
the  present  design  can  result  in  relatively  long  forebodies 
and  relatively  short  bases  with  steep  variation  in  local 
diameter,  a small  error  in  location  of  xBa  oould  result 
in  a large  error  in  db  if  x.a  was  measured  from  the  nose. 
Hence,  much  more  accurate  results  are  obtained  in  the  pres- 
ent case  by  measuring  xBa  from  the  beginning  of  the  base 
section.  Onae  xBa,  and  thereby  db,  is  determined  from 
the  figure,  the  potential  force  can  be  calculated  from 
equation  (12) . 


Viscous  Crossflow  Force . Crossflow  separation  usually 
begins  at'  the  base  of  the  body  when  an  angle  of  attack  of 
approximately  6°  is  reached.  With  increasing  angle  of  attack, 
the  location  of  this  crossflow  separation  occurs  further 
upstream,  so  that  eventually  most  of  the  body  is  in  separated 
crossflow.  A previous  correlation  of  available  data  (ref. 

8)  has  shown  that  the  location  of  the  beginning  of  cross- 
flow  separation  can  be  expressed  as 

■ 1 + (14) 
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The  value  of  the  viscous  force  term  can  be  written  as 


sin*a 

K 


(15) 


where  S„q  is  the  planform  area  of  that  portion  of  the  body 
in  separated  crossflow.  The  crossflow-drag  coefficient  is 
a function  of  the  local  crossflow  Reynolds  number.  There- 
fore, for  sections  of  changing  diameter  a different  value 
for  Cd  is  required  for  every  section.  Values  of  Cd 
c ac 

are  obtained  from  figure  18,  taken  from  reference  9.  These 
values  of  Cd  are  average  values  over  the  entire  separated 
c 

region.  This  is  all  that  is  necessary  for  the  computation 
of  normal  force.  However,  when  computing  pitching  moment, 
the  axial  distribution  of  crossflow  drag  is  required.  This 
point  is  discussed  in  more  detail  in  the  next  section. 

This  method  for  calculating  the  viscous  crossflow  force 
is  contained  in  equations  (14)  and  (15)  plus  figure  18.  It 
is  a well-established  and  documented  procedure  which  has 
been  adapted  to  bodies  of  the  present  type. 

A description  of  method  use  will  now  be  given,  followed 
by  a numerical  example. 

Use  of  Method 


The  parameters  required  to  determine  the  body-alone 
normal-force  coefficient  are 

Slenderness  ratio  of  forebody,  lp/d 

Local  body  radius  distribution,  r(x) 

The  procedure  to  be  followed  is: 

(1)  Given  the  base  and  1 «/d,  determine  xa.  from 

figure  17.  r aa 

(2)  From  xsa,  determine  du,  the  diameter  of  the 
base  at  xaa.  This  is  the  effective  base  of 
the  body. 
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(3)  Uae  djj  in  equation  (13)  to  calculate  CN  . 

a 

(4)  Find  C„  from  equation  (12) . 

Npot 

(5)  For  each  angle  of  attack  at  which  the  value  of 

CN  is  required,  determine  x>c  from  equation 

(14?. 

(6)  Knowing  the  location  of  crossflow  separation, 

determine  . 

pc 

(7)  Determine  Cd  from  figure  18  for  each  local 

Qc 

section  behind  xu_. 

Hu 

(8)  Calculate  CM  from  equation  (15). 

Nnl 

(9)  Add  the  results  from  steps  (4)  and  (8)  to  deter- 

mine total  CM  . 

nb 

Numerical  Example.  This  example  compares  the  method  to 
data  from'  the  present  test.  The  configuration  chosen  is 
N2C2B2 , at  a free  stream  Reynolds  number  of  35x10®.  The 
values  of  the  required  parameters  are 

Forebody  slenderness  ratios  1^/d  ■ 5 

Nose  sections  One  caliber  ellipsoid 

Base  sections  Two  caliber  circular  cone 

Body  diameters  7.0  in. 

A sketch  of  this  configuration  is  given  in  figure  19. 

(1)  From  figure  17,  xaa/lb  ■ 0.455. 

(2)  The  equation  for  the  diameter  of  a cone  two 

calibers  in  length  is 

dx  - d(l  - §) 

where  d is  the  diameter  of  the  base  of  the 
cone.  Thus  d^  ■ 3.815  in.  and  d^/d  ■*  0.545. 

(3)  From  equation  (13) , CN  ■ 0.594  par  radian. 

■*  Al 
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FIGURE  19.  BODY  CONFIGURATION  N2C2B2 
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(4) 

Use  C„  -■  C_,  a to  find  potential  force  as  a 

pot  Cl 

function  of.  angle  of  attack.. 

(5) 

Determine  xgc  from  equation  (14) 
of  attack  required. 

for  each  angle 

(6) 

Calculate  planform  area  from  each 
(5)  to  end  of  body. 

xBC  in  step 

(7) 

Find  the  appropriate  value  of  Cd 

for  each  local 

6 

section  at  each  crossflow  Reynolds 

number. 

In  the  present  example  Cd  * 0.290  throughout. 

Gc 

(8)  Use  equation  (15)  to  find  CM 

Nnl 

(9)  Add  result  from  step  (4)  to  result  from  step  (8) 

to  find  CM  . 

nB 


Steps  (4)  through  (9)  are  summarized  in  the  table  below. 


Step 

Step 

Step 

Step 

Step 

Step 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

a 

»,0/a 

pa 

cd 

a 

Cnb 

2 

. 

0.0207 

- 

- 

- 

0.0207 

- 

4 

0.0414 

- 

- 

- 

- 

0.0414 

6 

0.0621 

3.500 

0.8897 

0.290 

0,0106 

0.0727 

8 

0.0827 

2.250 

1.3537 

0.290 

0.0285 

0.1112 

10 

0.1031 

1.833 

1 . 5250 

0.290 

0.0499 

0.1530 

12 

0.1235 

1.625 

1.6234 

0.290 

0.0762 

0.1997 

14 

0.1437 

1.500 

1.6924 

0.290 

0.1075 

0.2512 

< 

16 

0.1637 

1.417 

1.7466 

0.290 

0.1440 

0.3077 

• ' 1 

This  is  compared  to  the  experiment  in  figure  20.  The  agree- 
ment is  quite  good  throughout  the  angle  of  attack  range. 
Further  comparisons  between  estimates  and  data  are  shown  in 
figures  21  and  22. 
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FIGURE  20.  COMPARISON  BETWEEN  PREDICTED  AND  EXPERIMENTAL 
BODY-ALONE  NORMAL  FORCE  COEFFICIENT. 
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FIGURE  21.  COMPARISON  BETWEEN  PREDICTED  AND  EXPERIMENTAL 
BODY-ALONE  NORMAL  FORCE  COEFFICIENT. 
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FIGURE  22.  COMPARISON  BETWEEN  PREDICTED  AND  EXPERIMENTAL 
BODY-ALONE  NORMAL  FORCE  COEFFICIENT. 
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V . 3 BODY-ALONE  PITCHING-MOMENT  COEFFICIENT 

This  method  permits  estimation  of  the  pitching-moment 
coefficient  on  isolated  bodies  at  angles  of  attack  up  to 
15°.  The  method  consists  of  a linear  term  derived  from 
potential-flow  theory  and  a nonlinear  term  which  is  the  con- 
tribution from  viscous  crossflow.  Included  in  the  nonlinear 
term  is  an  estimation  of  the  axial  distribution  of  the  cross- 
flow normal  force.  Comparison  with  data  from  the  present 
experiment  shows  reasonable  agreement  for  both  terms. 

An  adequate  estimation  of  the  pitching  moment  on  a 
slender  body  of  revolution  has  been,  in  the  past,  a very 
elusive  objective.  It  requires  not  only  a suitable  estimate 
of  the  normal  force,  but  also  a knowledge  of  the  axial  dis- 
tribution of  that  normal  force.  As  developed  in  the  pre- 
vious section,  the  normal  force  consists  of  a linear  term 
due  to  potential  lift  and  a nonlinear  viscous  term.  Thus, 
the  distribution  of  two  forces  are  actually  required.  How- 
ever, in  addition  this  means  that  the  pitching  moment  also 
can  be  treated  as  the  sum  of  a linear  term  (potential  flow) 
and  a nonlinear  term  (viscous  flow).  Thus,  the  objective 
of  this  section  is  to  develop  an  expression  for  pitching- 
moment  coefficient  of  the  form 


a,a  + a,a‘ 


(16) 


Applying  slander  body  theory,  Allen  and  Perkins  (ref. 
4)  found  the  following  expression  for  the  pitching  moment j 


It 


'Mr 


1 

*7 


dS 

9x 


(x-xm)dx( 


sin  2a  cos 


I) 


(17) 


'B 


Sb(1B  " V 


u 


sir,  2a  cos 


(18) 


This  analysis,  applied  to  bodies  that  are  closed  at  the 
base  ( S)-,  ■0),  as  ware  those  here,  yields  the  result  that 
the  pitching  moment  is  only  a function  of  the  size  (i.e., 
volume)  of  the  body.  But,  as  discussed  in  the  previous 
section  on  normal-force  coefficient,  there  is  an  effective 
base  due  to  the  separation  or  thickening  of  the  axial-flow 
boundary  layer.  The  axial  location  of  the  affective  base, 
xga,  marks  the  downstream  extant  of  the  potential  flow. 
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Thus,  it  seems  more  realistic  to  use  this  position,  rather 
than  the  end  of  the  body,  as  the  upper  limit  of  integration 
in  equation  (17) . 

The  usual  practice  of  estimating  the  contribution  to 
the  pitching  moment  due  to  viscous  crossflow  (e.g.,  see  ref. 
7 or  10)  is  to  assume  a relation  of  the  form 


- 2 


'M 


hi 


2naz 

V" 


B 


sc 


- x)dx 


(19) 


In  this  expression 


assumes  the  value  corresponding  to 


the  local  crossflow  Reynolds  number.  For  a cylindrical 
body  then,  Cd  is  usually  assumed  constant  with  x and  equa- 
c 

tion  (19)  gives  a viscous  force  distribution  that  is  uniform 
in  x.  The  result  of  this  assumption  usually  is  that  the 
nonlinear  pitching  moment  is  under-predicted  (ref.  7).  For 
this  reason  a brief  analysis  of  the  expected  form  of  the 
viscous  force  axial  distribution  is  helpful. 

If  the  separated  crossflow  boundary  layer  is  represented 
by  a set  of  point  vortices,  then  the  vortex  impulse  theorem 
for  moderate  angles  of  attach  gives  the  following  expression 
for  the  loading  on  an  infinitesimal  length  of  the  body  In 
the  presence  of  vortices  (ref.  11) 


Nv 

dN  * Real  l d(To)  . 

j-1  3 


(20) 


whore  T is  the  vortex  strength,  a is  the  complex  distance 
between  a vortex  and  its  image  inside  the  cylinder,  and  N 
is  the  number  of  vortices.  For  moderate  angles  of  attack 
the  lee-side  vortex  wake  normally  can  be  represented  by  two 
vortices  which  reduces  equation  (20)  to 


dN  « d(Tbv)  - bydr  + 


rdb 


(21) 


where  bv  is  the  distance  between  the  vortices  or  the  span 
of  the  vortex  wake.  Thus,  the  axial  distribution  of  the 
vortex-induced  force  becomes 
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dN 

3Z 


dr 


d7  + 


db 

v 

■ar 


(22) 


Measurements  of  both  by  and  T have  been  correlated  In 
reference  8 by  Mendenhall  and  Nielsen , who  show  that,  to 
a first  approximation,  both  terms  are  linear  in  x.  These 
results  are  reproduced  in  figures  23  and  24. 

Since  dbv/dx  is  small  it  can  bo  ignored  in  equation 
(22)  so  that 


dN 


ax 


oc 


dr 

ax 


And,  since  r * x,  then  the  normal  force  will  be  dis- 
tributed in  an  approximately  linear  fashion  along  the  body. 
Therefore,  for  purposes  of  estimating  the  pitching  moment 
due  to  the  viscous  normal  force,  the  crossflow-drag  coeffi- 
cient can  be  represented  by  the  relation 

V"  - (sr^K 

This  relation  distributes  the  crossflow-drag  force  linearly 
from  the  position  of  crossflow  separation  xso  to  the  end 
of  the  body.  The  crossflow-drag  coefficient,  Cd  , is  the 

c 

value  based  on  local  crossflow  Reynolds  number  and  is  given 
in  figure  18. 

Description  of  Method 

Pitohing-Momsnt-Curve  S lope . The  present  method  in- 
volves  an  equation  of  the  form  given  in  equation  (16)  and 
repeated  here  for  convenience. 

c„b  . ...  + «..* 

The  one  boundary  condition  available  is 

CM  (a  - 0)  - 0 (24) 
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FIGURE  23.  BODY-VORTEX  CHARACTERISTICS  - VORTEX 
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Differentiating  equation  (16)  with  respect  to  x yields  the 
result  that 


i.e.,  a 2 is  the  pitching-moment-curve  slope.  In  the  pres- 
ent method,  the  following  expression  is  used  to  estimate 


K_(K, 

n\ 


*i> 


sa 


TF 


dS 

3? 


(x_  - x)dx 
m 


(26) 


The  term  (K2  - K l ) is  the  apparent  mass  factor  of  Munk 
(ref.  5)  for  ellipsoid-shape  bodies  and  is  given  in  figure 
25.  The  term  K is  an  empirical  factor  derived  from  the 
data  and  appliedmto  the  apparent  mass  term.  This  term  is 
shown  in  figure  26.  The  upper  limit  of  integration  is  the 
axial  loaation  of  the  effective  base.  This  quantity,  along 
with  a more  aomplete  discussion  of  the  procedure  used  to 
obtain  its  value  can  be  found  in  Section  V.2. 

Nonlinear  Pitching  Moment.  It  is  assumed  that  the 
viscous  dross^low  exerts  a force  on  the  body  from  xiC,  the 
axial  location  of  crossflow  separation,  to  the  end  or  the 
body,  and  that,  along  this  length,  the  viscous  force  varies 
linearly  with  x.  With  these  assumptions  the  following  term 
represents  the  pitching  moment  due  to  viscous  crossflow. 


The  crossflow  drag  coefficient  is  the  value  corresponding  to 
the  local  crossflow  Reynolds  number  at  each  section  of  the 
body  between  x8_  and  xQ.  Curves  for  C.  are  given  in 

c 

figure  18.  The  term  n is  the  ratio  of  the  drag  on  a 
finite  cylinder  to  the  drag  on  on  infinite  cylinder  and  is 
applied  as  a correction  factor  for  finite  length.  Values 
for  n are  obtained  from  data  given  by  Goldstein  in  refer- 
ence 12  and  are  presented  in  figure  27.  The  lower  limit  of 
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FIGURE  27.  RATIO  OF  CROSSFLOW-DRAG  COEFFICIENT  FOR  A FINITE 
LENGTH  CYLINDER  TO  THAT  FOR  AN  INFINITE  LENGTH 
CYLINDER  (FROM  REF.  12). 


integration,  xao,  is  given  by  the  following  expression,  based 
on  a correlation  of  experimental  data  by  Mendenhall  and 
Nielsen  in  reference  8. 


x 


sc 


~r 


(28) 


where  a is  in  degrees. 

Use  of  Method 

The  following  parameters  are  required  in  the  applica- 
tion of  this  method. 

Body  slenderness  ratio,  l0/d 

Forebody  slenderness  ratio,  lp/d 

Distribution  of  cross-sectional  area  along  body,  S(x) 
Position  of  the  moment  reference  center, 

m 

Local  crossflow  Reynolds  number,  Re  (x) 

n 

Axial  distribution  of  radius,  r(x) 

Length  of  body,  1B(«  x0) 

(1)  Determine  the  apparent  mass  factor  (Ka  - Ki),  appro- 

priate for  the  particular  configuration,  from  figure  23. 

(2)  Determine  K from  figure  26. 

m 

(3)  Find  x,,,  from  figure  17. 

BA 

(4)  Use  equation  (26)  to  determine  the  pitching-moment- 

curve  slope. 

(3)  Determine  n from  figure  27. 

(6)  For  each  angle  of  attack  required,  use  equation  (28) 

to  find  the  value  of  x.„. 

sc 

(7)  Find  the  local  value  of  Cd  fer  every  section  from 

figure  18.  °c 

(8)  Use  equation  (27)  to  determine  the  nonlinear  pitching 

moment . 
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Numerical  Example.  This  example  compares  the  method  with 
data  from  the  present  test.  The  configuration  chosen  is 
N2C2B2  at  a free  stream  Reynolds  number  of  16.2x10s , based 
on  body  length.  The  planform  of  this  configuration  is 
given  in  figure  19.  The  values  of  the  required  parameters 
are 


Body  slenderness  ratio t lg/d  » 7. 

Forebody  slenderness  ratio:  lF/d  ■ 5. 

Distribution  of  cross-sectional  area  along  body: 

This  body  can  be  described  analytically.  The  following 
equations  desaribe  the  cross-sectional  area. 

Nose:  S (x)  « J [d2  - (x  - d)  2] 

ltd* 

Center  seation:  S(x)  ■ -j- 


0 < x < d 
d < x £.  5d 


Base: 


S(x) 


(x  - 5d) 
2d 


5d  < x S.  7d 


Body  diameter:  d ■ 7 in. 


Position  of  moment 

reference  center:  x_  ■ 30.5  in. 


Axial  distribution  of  radius 

Nose:  r{x)  - ^ — — -g— IL.j 

Center  section:  r(x)  ■ ^ 

Base:  r(x)  - | Jl  - (* 

Length  of  bodyi  xg  ■ 49  in. 


1/2 


0 L x £ d 
d < x £ 5d 
5d  < x £ 7d 
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(1)  Prom  figure  25,  (Ka  - Ki)  = 0.898 

(2)  Prom  figure  26,  Km  * 0.940 

(3)  Using  figure  17,  x_a/lv  ■ 0.455.  Since  l^/d  * 2 for 

base  B2 , xBa/d  - *7910,  measured  from  the  start  of 
the  base  section.  Measured  from  the  nose, 
x_a/d  - 5.91  ■+■  xga  ■ 3.4475  ft.  This  is  the  upper 
limit  of  integration  in  equation  (26)  . 

(4)  Substituting  the  necessary  quantities  into  equation 

(26)  and  carrying  out  the  integration 


'M. 


- 8.043  rad 


_ i 


(5) 

(6) 


From  figure  27,  n - 0.635 

Find  x._/d  at  7°,  10°,  13°,  from  equation  (28) 


(7) 

(8) 


ot 

7 

10 

13 


xec'd 


2.667 

1.U3  3 

1.556 


The  crossflow-drag  coefficient,  from  figure  18  is  0.290 
for  all  sections. 

The  results  of  the  integration  of  equation  (27)  are 

CMnl 


a 

7 

10 

13 


-0.4367 

-0.0091 

+0.0963 


Thus,  the  final  result  for  the  pitching  moment  is  shown  in 
figure  28.  The  agreement  is  quite  good,  especially  at  the 
lower  angles  of  attack.  This  indicates  that  the  estimation 
of  th£  pitahing-moment-curve  slope  is  very  accurate.  The 
present  method  does  underestimate  the  nonlinear  portion  of 
the  pitching  moment  somewhat,  but  for  initial  design  purposes 
the  agreement  is  considered  adequate.  Possibilities  for  the 
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FIGURE  23.  COMPARISON  BETWEEN  PREDICTED  AND  EXPERIMENTAL  BODY' 
ALONE  PITCHING-MOMENT  COEFFICIENT. 
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causes  of  the  discrepancy  in  Cj.  include  the  position  of 

nl 

crossflow  separation,  x.c,  and  the  form  of  the  viscous-force 
distribution  along  the  body  axis.  Additional  comparisons 
are  given  in  figures  29  and  30.  Agreement  is  quite  good  for 
both  cases. 


V . 4 BODY-ALONE  AXIAL-FORCE  COEFFICIENT 


This  method  permits  estimation  of  the  body-alone  axial- 
force  coefficient  on  isolated  bodies  in  incompressible  flow 
at  angles  of  attack  up  to  15°.  This  procedure  consists  of 
first  determining  the  zero  angle  axial  force  and  then  the 
slope  of  its  curve.  Comparison  with  data  from  the  present 
test  shows  very  good  agreement  for  all  angles. 

The  axial-force  coefficient  on  isolated  bodies  of  rev- 
olution is  typified  by  a finite  value  at  zero  angle  of 
attack  (zero-angle  drag) . For  bodies  used  in  the  present 
test  the  behavior  as  the  angle  of  attack  increased  from 
zero  was  approximately  linear  with  |a|.  Hence,  a repre- 
sentative equation  is  one  of  the  form 


(29) 


For  incompressible  flow,  the  sources  of  zero-angle 
axial  force  are  base  drag  and  skin  friction;  i.e., 


a 


o 


+ C_ 


(30) 


It  is  common  practice  (e.g. , refs.  10  and  13)  to  treat  these 
two  sources  independently,  Drag  due  to  skin  friction  is 
based  on  the  well-documented  values  for  turbulent  skin- 
friction  coefficient  on  a flat  plate  (e.g.,  ref.  13), 
modifying  for  thickness  effects  of  a body  of  revolution, 
and  applying  it  to  the  actual  surface  area  of  the  body 
under  test.  Based  on  experimental  data,  plus  the  best 
available  information  for  calculating  skin  friction  (ref. 

13) , equation  (30)  was  used  to  identify  the  base  drag  in 
the  present  test. 


Examination  of  the  data  revealed  a strong  dependence 
of  C on  base  length.  In  addition,  the  base  drag  for  a 


particular  configuration  displayed  a slight  dependence  on 
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the  free  stream  Reynolds  number,  being  slightly  greater  at 
the  larger  Reynolds  number.  There  was  also  a dependence  on 
the  length  of  the  body  that  can  be  related  to  skin-friction 
effects.  Finally,  negligible  effect  of  nose  type  was 
observed. 

Description  of  Method 

Zero-Angle  Axial  Force,  (C^  } . To  determine  the  value 

B 

of  the  axial-force  coefficient  at  zero  angle  of  attack  the 
contributions  from  the  base  and  skin  friction  are  summed. 

For  bodies  with  conical  base  sections  and  sharp  center 
section-base  section  junctions,  the  base  drag  is  given  in 
figure  31  as  a function  of  base  length.  The  curves  in  thiB 
figure  were  determined  as  described  above  from  equation  (30) 
using  zero-angle  axial  force  for  three  configurations  having 
identical  forebodies  and  differing  bases. 

Drag  due  to  skin  friction  on  bodies  of  revolution  is 
well  predicted  using  values  for  the  turbulent  skin-friction 
coefficient  on  a flat  plate.  This  value  is  applied  to  the 
actual  surface  area  exposed  to  potential  flow  (i.e.,  nose, 
downstream  to  xsa) , with  a teem  for  thickness  effeots 
included.  Suah  a relation  is  expressed  by  the  following 
equation. 


The  quantity  in  brackets  is  a geometrical  term  for  the 
effect  of  thickness  (ref.  13).  Normally,  the  length  chosen, 
Isa'  the  distance  back  to  x-_,  the  point  of  axial  flow 
separation.  Similarly,  S is  tne  surface  or  wetted  area 
from  the  nose  to  xBa.  Tne  value  of  Cf,  the  skin-friction 
coefficient  based  on  wetted  area,  is  given  by  the  Schoenherr 
equation  (ref.  13  or  14),  equation  (32),  and  plotted  in 
figure  32  for  convenience.  The  equation  is  based  on  a 
similarity  rule  derived  by  von  Karman  (sea  ref.  13). 


log(Re-aCf)  - 


0.242 

■*? 


(32) 


In  this  relation,  Cf  is  the  turbulent  skin-friction  coeffi- 
cient based  on  exposed  surface  area  and  Reaa  is  the 
Reynolds  number  based  on  plate  length.  Thus,  when  deter- 
mining tho  value  of  Cf  from  figure  32,  the  value  of 
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Reynolds  number  based  on  the  length  of  the  body  exposed  to 
potential  flow  dsa)  should  be  used. 

Variation  of  C*  with  Angle  of  Attack.  Analysis  of 

ab 

the  data  from  the  present  teat  showed  that,  to  a very  close 
approximation,  Ca  varied  linearly  with  angle  of  attack. 

B 

The  value  of  the  axial-force-curve  slope  is  essentially 
independent  of  nose  shape  and  Reynolds  number  and  is  basi- 
cally a function  of  the  length  of  the  center  section  plus 
the  base  section.  The  variation  of  the  axial-force-curve 
slope,  aj,  with  center  section  and  base  section  length  is 
shown  in  figure  33  which  is  based  on  the  experimental  data. 

Use  of  Method 


This  section  demonstrates  how  the  method  is  to  be  used 
and  makes  comparisons  with  data.  The  following  parameters 
are  required  for  the  application  of  this  method. 

Base  length,  in  calibers,  1^/d 

Length  of  body  exposed  to  potential  flow,  in 
calibers,  1/d 

S a 

Surface  area  of  portion  of  body  exposed  to  potential 
flow,  Sw. 

Reynolds  number  based  on  the  length  l8a, 

Center  section  length,  in  calibers,  l_/d. 

v 


General  Description. 

1.  Determine  base  drag,  C 


, from  figure  31. 


2.  Find  skin-friction  coefficient  from  figure  32  or 

equation  (32). 

3.  Calculate  skin-friction  drag,  C , from  equation 

(31). 


4.  Add  results  from  steps  (1)  and  (3)  to  find  zero- 

angle  axial  force,  (C.  ) . 

mB  o 

5.  Find  axial-force-curve  slope,  air  from  figure  33. 


| 

1 
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6.  Use  equation  (29)  to  find  axial-force  coefficient 
as  a function  of  angle  of  attack. 

Numerical  Example.  This  example  compares  the  method  to 
data  from  thepresenttest . The  configuration  chosen  is 
N3C2B2,  at  a free-stream  Reynolds  number  of  40x10s  based 
on  body  length.  The  values  of  the  required  parameters 
are  t 


• Base  length,  1^/d  ■ 2.0. 

• Length  of  body  exposed  to  potential  flow,  using 

figure  17  (Section  V.2),  l_,/d  « 6.860. 

BA 

• Surface  area  of  body  exposed  to  potential  flow, 

S,,  - 6.743  ft. 
w 

• Reynolds  number  based  on  1__,  Re„.  *•  34.18x10s , 

Bel  Ba 


• Center  section  length,  lc/d 


1.  From  figure  31,  C 


6.0. 

0.0282. 


2.  Figure  32  gives  C ^ * 0.00242. 

3.  Using  equation  ("31)  , Cn  ■ 0.0661. 

- 0.0943. 


4.  Adding  C and  C 


'f 

<CA> 


B o 


5.  From  figure  33,  a2  ■ 0.1089  per  radian. 

6.  Thus,  C . » 0.0943  + 0.1089  a 

aB 


This  estimate  is  compared  to  the  experimental  values  in 
figure  34.  The  agreement  is  excellent  throughout  the  angle 
of  attack  range.  Further  comparisons  between  the  experi- 
ment and  estimates  from  the  present  method  are  shown  in 
figures  35  and  36.  Agreement  is  excellent  in  both  cases. 
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VI.  ISOLATED  TAILS 


VI. 1 GENERAL  CONSIDERATIONS 

When  a low  aspect  ratio  tail  is  placed  in  an  airatream, 
attached  and  separated  flows  are  induced,  in  much  the  same 
way  as  on  a body  alone.  At  low  angles  of  attack,  the  flow 
is  attached  over  the  whole  tail.  At  angleB  greater  than 
about  6 degrees,  however,  the  boundary  layers  at  the  lead- 
ing  end  side  edges  separate  and  roll  up  into  large  vortices 
on  the  upper  surface,  Big.  37.  These  vortices  produce 
suction  which  increases  the  total  normal  force  on  the  wing. 
This  normal  forae  increases  nonlinearly  with  angle  of 
attack.  The  magnitude  of  the  nonlinear  force  has  been 
found  (Ref.  1)  to  be  dependent  upon  several  of  the  wing 
geometrical  parameters.  In  particular,  the  force  increases 
as  the  leading  edge  becomes  sharper.  For  tails  typical  of 
submersible  vehicles,  the  relatively  blunt  leading  edges 
used  moan  that  additional  control  power  could  be  obtained 
by  sharpening  the  leading  edge. 


VI. 2 TAIL-ALONE  NORMAL-FORCE  COEFFICIENT 

This  method  permits  calculation  of  the  normal-force 
coefficient  of  low-aspect-ratio  isolated  tails  with  roundod 
leading  edges  and  blunt  trailing  edges.  This  method  is 
valid  for  incompressible  flow  and  for  angles  of  attack 
up  to  30°.  The  method  is  divided  into  two  steps.  First, 
the  normal-force  curve  slops  at  zero  angle  is  determined  as 
a function  of  aspect-ratio  and  trailing-edge  shape.  Second, 
a nonlinear  term  is  determined  as  a function  of  aspect 
ratio,  taper  ratio,  and  Reynolds  number  to  account  for 
viscous  effects. 

The  normal  force  on  a low-aspect-ratio  lifting  surface 
(wing  or  tail)  can  be  thought  of  as  consisting  of  two  con- 
tributions. The  first  is  well-approximated  by  potential-flow 
theory  and  is  linear.  The  second  is  due  to  the  action  of 
viscosity  and  is  highly  nonlinsar.  At  angles  of  attack 
greater  than  a few  degrees,  tail  normal  force  o&n  become 
dominated  by  the  nonlinear  viscous  effects.  Potential-flow 
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theory  does  not  take  the  nonlinear  phenomena  into  account. 

It  would  be  desirable  to  be  able  to  calculate  both  of  these 
effects  exactly.  However,  the  current  state  of  nonlinear 
theory  is  such  that  it  is  not  possible  to  do  this,  in 
general,  with  any  accuracy.  Fortunately,  for  preliminary 
design  purposes,  an  exact  calculation  is  unnecessary.  The 
following  semi  empirical  method  was  formulated  to  combine 
the  two  effects  in  a convenient  manner  for  use  in  prelim- 
inary design. 

The  method  is  similar  to  those  based  on  the  crossflow 
drag  analogy  such  as  described  by  Flax  and  Lawrence  (Ref. 

15)  . In  these  methods  the  normal-force  coefficient  is 
assumed  to  be  of  the  form 

Cxi“aja  + aj,a2  (33) 

WT 

Taking  the  derivative  of  this  equation  we  see  that  at  is 
equal  to  the  zero  angle  normal-force  curve  slope,  i.e., 


The  constant  a2  is  chosen  such  that  the  full  expression  pro- 
vides the  best  overall  comparison  to  the  experimental  data. 

Normal-Force  Curve  Slope 

Correlation  of  the  experimental  data  showed  that,  for 
aspect  ratios  less  than  or  equal  to  one,  the  normal-force 
curve  slope  is  very  well  predicted  by  slender  wing  theory; 
i.e., 


CN  - J AR  for  ffi  < 1.0  (35) 


For  higher  aspect  ratios  the  experimental  data  fell  below  the 
slender  body  theory  prediction.  This  trend  is  well  known  and 
is  documented  in  reference  15.  More  refined  potential  flow 
theories  are  available  which  predict  this  trend  and  give 
much  better  agreement  with  experiment  at  higher  aspect  ratios 
for  wings  with  sharp  trailing  edges.  A comparison  between 
data  from  the  current  work  and  from  previous  tests  is  shown 
in  figure  38.  It  may  be  seen  that  the  normal-forca-curve 
slopes  derived  from  the  present  data  are  generally  greater 
than  those  found  in  other  experiments.  This  is  believed  to 
be  an  effect  of  trailing-edge  shape  as  shown  by  Hoerner 
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FIGURE  38.  MEASURED  VALUES  OF  LIFT-CURVE  SLOPE. 
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(Ref.  16).  A wing  section  with  a blunt  trailing  edge,  as  in 
the  present  experiment,  has  a greater  normal- force-curve 
slope  than  the  same  section  with  a sharp  trailing  edge,  as 
seen  in  figure  39  (Ref.  16).  Hence,  a suitable  prediction 
method  for  should  include  trailing-edge  shape  as  a 

parameter.  a 

Examination  of  the  data  for  the  effect  of  taper  ratio 
shows  that,  for  blunt  trailing  edges  (the  present  data) , 
taper  ratio  has  a negligible  effect.  For  sharp  trailing 
edges  however,  there  is  a slight  effect,  with  rectangular 
wings  having  slightly  higher  normal-force-curve  slopes  than 
delta  wings. 

The  effect  of  Reynolds  number  over  the  range  tested 
(4x10®  to  8.5xlo6  per  foot)  had  an  insignificant  influence 
on  C„  . This  was  confirmed  by  a study  by  Jones  (Ref.  17) 

Na 

who  found  only  a very  small  increase  in  C with  Reynolds 
number . Na 

It  is  clear  then,  that  any  method  for  predicting  normal- 
force-curve  slope  should  include  aspect  ratio,  shape  of  the 
trailing  edge,  and,  for  wings  with  a sharp  trailing  edge, 
the  taper  ratio.  These  effects  are  summarized  in  figures 
40 (a)  and  (b)  which  constitute  the  recommended  curves  for 
determination  of  the  normal-force-curve  slope.  Figure  40(a) 
is  for  blunt  trailing  edges  and  is  a best  fit  curve  to  the 
present  data.  For  aspect  ratios  less  than  or  equal  to  one 
it  coincides  with  the  value  obtained  from  slender-wing- 
theory.  The  theoretical  curves  of  figure  40(b),  recommended 
for  sharp  trailing  edges,  are  due  to  Lawrence  (Ref.  IB). 

This  theory  provides  adequate  approximations  to  lifting- 
surface  theory  for  plariforms  with  straight  trailing  edges 
of  aspect  ratio  less  than  three. 

Nonlinear  Effects 


The  understanding  and  prediction  of  nonlinear  effects 
on  low-aspect  ratio  wings  is  not  as  advanced  as  the  linear 
portion.  Much  early  work  was  done  trying  to  develop  an 
expression  in  terms  of  a cross-flow-drag  coefficient.  This 
theory  relates  the  nonlinear  effects  to  the  drag  associated 
with  the  two-dimensional  viscous  flow  around  an  infinitely 
long  flat  plate  in  a flow  having  a velocity  equal  to  the 
component  of  the  f ree-stream  velocity  normal  to  the  wing . 
This  concept,  which  has  been  employed  here,  requires  the 
introduction  of  a suitable  empirical  drag  coefficient. 
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Values  for  this  coefficient  were  selected  which  gave  the 
best  overall  fit  to  the  experimental  data. 

Correlation  of  the  present  data  resulted  in  the  formula 
tion  of  the  curves  presented  in  figures  41(a),  (b)  and  (c) 
which  serve  as  the  recommended  values  for  a 2 in  equation 
(33) . The  data  revealed  that  a2  depends  strongly  on 
aspect  ratio,  as  well  as  on  the  taper  ratio  and,  to  a lesser 
extent,  on  the  Reynolds  number.  Bartlett  and  Vidal  (Ref. 

19)  demonstrated  further  that  the  value  of  this  coefficient 
varies  with  the  shape  of  the  leading  edge,  being  greater 
for  sharp  leading  edges  than  for  rounded  leading  edges. 

Since  the  tails  in  the  present  experimental  study  had 
rounded  leading  edges,  the  curves  in  figure  41  should  be 
used  for  wings  with  similar  leading  edges. 

Use  of  Method 


This  section  demonstrates  the  use  of  the  method  and 
makes  comparisons  with  experimental  data.  The  required 
physical  characteristics  of  the  tail  are  the  following: 

Aspect  ratio  (based  on  complete,  symmetric  planform) 

Taper  ratio 

Shape  of  trailing  edge 

Reynolds  number  based  on  root  chord,  ReCR 

In  addition,  the  leading  edge  should  be  rounded  and  the  flow 
should  be  incompressible. 

General  Description. 

(1)  Determine  the  value  of  CN  from  the  appropriate 

curve  of  figure  40.  a 

(2)  Determine  the  appropriate  value  of  a2  from 
figure  41. 

(3)  Using  these  two  quantities  in  equation  (33) , find 
CM  an  a function  of  angle  of  attack  (a  in 

WT 

radians) . 

Numerical  Example  I . This  example  compares  the  method 
with  one  'of  the  cases  in  the  present  experimental  study. 

The  case  selected  is  tail  T10. 
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FIGURE  41(0).  RECOMMENDED  CURVES  FOR  PREDICTION  OF  a, , 

COEFFICIENT  OF  NONLINEAR  TERM  OF  TAIL-ALONE 
NORMAL  FORCE  COEFFICIENT;  X ■ 0.0. 
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Aspect  ratio  = 0.5 

Taper  ratio  » 0.5 

Shape  of  trailing  edges  Blunt 

Reynolds  number  ■ 4.35x10*,  based  on  the  root  chord. 

_ 2 

(1)  From  figure  40(a),  for  JR  « 0.5,  CN  ■ 0.7854  rad 

Nct 

(2)  From  figure  41  (bl,  for  JR  * 0.5,  and  X *•  0.5, 

aa  ■ 1.641  rad"2 

(3)  Applying  equation  (33)  we  calculate  the  following 

values. 


2° 

C»T 

4 

0.0628 

8 

0.1416 

12 

0.2364 

16 

0.3472 

20 

0.4740 

24 

0.6168 

28 

0.7756 

The  comparison  between  prediction  and  data  are  presented  in 
figure  42.  As  can  be  seen,  agreement  is  excellent.  This 
might  have  been  expected  since  the  data  themselves  formed 
the  basis  for  the  method.  The  comparison  is,  however,  a 
necessary  step  in  checking  the  self-consistenoy  of  the 
method. 

Numerical  Example  II.  This  example  compares  the  present 
method  with  data  of  Bartlett  and  Vidal  (Ref.  19)  on  round 
edge  delta  wings. 

Aspect  ratio  ■ 1.5 

Taper  ratio  ■ 0.0 


8hape  of  trailing  edges  Sharp 
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Reynolds  number  ■ 3xl06 , based  on  root  chord 

(1)  Prom  figure  40(b),  for  £R  - 1.5  and  X ■ 0.0, 

CM  ■ 1.74  rad-1 

Na 

(2)  Prom  figure  41(c),  for  Jft  - 1.5  and  X » 0.0, 
a2  - -0.094  rad“s 

(3)  Applying  equation  (33)  we  calculate  the  following 
values. 


£ 

4 

0.1207 

8 

0.2406 

12 

0.3595 

16 

0.4774 

20 

0.5945 

24 

0.7106 

28 

0.8259 

This  table  is  plotted  in  figure  43  along  with  the  experi- 
mental data.  The  prediction  is  in  reasonably  good  overall 
agreement  with  the  data.  Below  20s  it  is  very  good.  In 
particular,  the  prediction  of  CN  is  excellent. 

a 

Figures  44  and  45  present  similar  comparisons  of  normal- 
force  coefficients  for  delta  wings  of  aspect  ratio  2.0  and 
1.33,  respectively,  the  data  being  taken  from  references 
19  and  20. 


VI. 3 TAIL-ALONE  CHORDWISE  CENTER  OF  PRESSURE 

This  method  permits  estimation  of  the  chordwise  loca- 
tion of  center  of  pressure  for  low-aspect-ratio  isolated 
tails  with  rounded  leading  edges.  The  method  is  valid  for 
angles  of  attack  up  to  30°  and  for  inoompreasible  flow. 

It  is  derived  through  correlation  of  the  data  from  the 
present  test, 
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FIGURE  43.  COMPARISON  BETWEEN  PREDICTED  AND  EXPERIMENTAL 
TAIL-ALONE  NORMAL  FORCE  COEFFICIENT. 
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PiaURE  44.  COMPARISON  BETWEEN  PREDICTED  AND  EXPERIMENTAL 
TAIL-ALONE  NORMAL  FORCE  COEFFICIENT. 


NCSC  TM-238-78 


NCSC  TM-238-78 


Chordwise  center  of  pressure  locations,  Xg«,  were  calcu- 
lated from  the  normal-force  and  hinge-moment  data.  The  loca- 
tions are  presented  in  terms  of  distance  from  the  root  chord 
leading  edge  and  are  nondimensionalized  by  the  root  chord 
CR.  For  angles  of  attack  greater  than  a few  degrees  the 
data  varied  smoothly  to  the  maximum  angle  of  attack_tested. 
Estimates  of  Xp  at  zero  degrees  angle  of  attack  (XT)Q 
were  made  by  determining  the  slopes  of  the  normal -force  and 
hinge-moment  coefficient  curves  at  zero  degrees.  These 
curves  formed  the  basis  for  the  method. 

The  zero-angle  centers  of  pressure  are  compared  to 
various  linear  theories  and  other  data  in  figure  46. 

Centers  of  pressure  from  the  present  investigation  are  seen 
to  be  slightly  downstream  of  both  theoretical  values  and 
measured  values  from  other  experiments . This  is  believed 
caused  by  the  blunt  trailing  edge  on  the  fins  used  in  the 
present  investigation.  As  discussed  in  the  previous  section, 
a wing  with  a blunt  trailing  edge  has  a greater  value  of 
the  normal-force  curve  slope  than  the  same  wing  profile  with 
a sharp  trailing  edge.  This  implies  additional  loading  over 
the  aft  portion  of  the  wing.  Hence,  a center  of  pressure 
that  is  further  downstream  than  a wing  with  a sharp  trail- 
ing edge  is  consistent  with  the  present  data. 

For  the  case  of  incompressible  flow  we  would  expect 
the  location  of  the  chordwise  center  of  pressure  to  be,  in 
general,  a function  of  taper  ratio,  aspect  ratio,  angle  of 
attack,  and  Reynolds  number;  i.e., 

XT  - 3^(1,*,  a,ReCR)  (36) 

Close  examination  of  the  data  revealed  that  XT  depended 
Btrongly  on  taper  ratio  and  angle  of  attack,  varied  slightly 
with  aspect  ratio,  and  wag  essentially  independent  of  the 
Reynolds  number.  These  findings  can  be  expressed  as 

XT  - Fj  U,a)Fa  (JR)  . (37) 

Curves  of  F,(A,a)  arc  shown  in  figures  47(a),  (b)  and  (c) 
for  the  three  aspect  ratios  tested.  The  differences  between 
the  curves  at  the  same  taper  ratio  shows  the  dependence  on 
aspect  ratio. 
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FIGURE  46. 


TAIL-ALONE  CHORDWISE  CENTER  OF  PRESSURE  AT 
ZERO  ANGLE. 
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FIGURE  47(a).  TAIL-ALONE  CHORDWISE  CENTER  OP  PRESSURE 
PREDICTION  CURVE > A R ■ 0.5, 
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Descr  i pt  ion  of  Method 

Linear  theory  predicts  no  variation  in  the  center  of 
pressure  location  with  angle  of  attack.  Fidler  and  Bateman 
(ref.  3)  found  a more  satisfactory  method  for  isolated 
tails  at  higher  Mach  numbers  to  be  one  based  on  the  corre- 
lation of  actual  data.  This  type  of  method  also  was  found 
to  bo  more  satisfactory  for  the  incompressible  flow  case. 
Thus,  the  method  consists  of  using  the  data  correlation 
curves  of  figure  47  discussed  earlier. 

Use  of  Method 

This  section  demonstrates  the  use  of  the  method  and 
makes  comparisons  with  experimental  data.  The  required 
physical  characteristics  of  the  tail  are  the  following: 

Taper  ratio 

Aspect  ratio  (based  on  complete  symmetric  planform) 

In  addition,  it  is  assumed  that  the  leading  edge  is  rounded 
and  the  flow  is  incompressible. 

General  Description . For  the  particular  combination 
of  taper  ratio  and  aspect  ratio  determine  the  variation  of 
xij.  with  angle  of  attack  from  the  appropriate  curve  of 
figure  47.  Use  linear  interpolation  for  aspect  ratio 
effects . 

Numerical  Example  J..  This  example  compares  the  method 
with  one  of  the  cases  "in  the  present  experimental  study. 
This  is  done  in  order  to  show  the  self-consistency  of  the 
method.  The  case  selected  is  tail  T10. 

X = 0.5 

PR  = 0. 5 

From  figure  47(a)  (aspect  ratio  0.5)  determine  Xm  as  a 
function  of  angle  of  attack  from  the  taper  ratio  0.5  curve. 
This  gives  the  following  values. 


0 0.3875 

5 0.3975 

10  0.4350 
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ft  ° 

X 

T/CR 

15 

0 

.4  650 

20 

0 

.4875 

25 

0 

. 5000 

30 

0 

. 5050 

The  comparison  between  prediction  and  data  is  shown  in 
figure  48.  As  can  be  seen,  agreement  is  excellent. 

Numerical  Example  11.  This  example  compares  the 
method  to  data  of  Bartlett  and  Vidal  (ref.  19)  on  round 
edge  delta  wings.  Also,  it  demonstrates  the  procedure  to 
be  applied  to  determine  the  effect  of  aspect  ratio. 

A = 0.0 

AR  = 1.5 

Determine  XT  as  a function  of  angle  of  attack  from  both 
figure  47(b)_,  (AR  = 1.0),  and  figure  47(c),  (AR  = 2.0). 
Then  find  Xm(a)  for  an  aspect  ratio  of  1.5  by  interpolat- 
ing between  these  values. 


0 

a_ 

xt/cr(ar  = 1,0) 

Xt/Cr(AR  = 2.0) 

xt/cr(ar  = 

0 

0.6375 

0.5875 

0.6125 

5 

0.6375 

0.5700 

0.6038 

10 

0.6275 

0.5650 

0.5963 

15 

0.6250 

0.5600 

0.5925 

20 

0.6225 

0.5600 

0 . 5913 

25 

0.6225 

0.5550 

0.5888 

30 

0.6225 

0.5425 

0.5825 

The  predicted  values  are  plotted  in  figure  49  with  the 
experimental  values.  Agreement  between  prediction  and 
experiment  is  seen  to  be  excellent  with  center  of  pressure 
location  predicted  within  2 percent  of  the  root  chord 
throughout  the  angle  of  attack  range. 
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Figure  50  shows  the  results  of  an  additional  comparison 
for  a delta  wing  of  aspect  ratio  2.0  and  the  data  from 
several  different  experiments  (refs.  19,  20  and  21).  The 
predicted  center  of  pressure  is  slightly  upstream  of  the 
experimental  values  but  still  is  within  2-1/2  percent  of 
the  root  chord  of  the  average  experimental  value . 
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o 

Experiment 

(Ref.  19) , 

ReCR  - 3«10‘ 

□ 

Experiment 

(Ref.  19) , 

ReCR  ■ 6«10‘ 

A 

Experiment 

(Ref.  20) 

o 

Experiment 

(Ref.  21) 

1 

Method 

I 

1 

i i i i 

a,  Degrees 


FIGURE  SO.  COMPARISON  BETWEEN  PREDICTED  AND  EXPERIMENTAL 
TAIL-ALONE  CHORDWI8E  CENTER  OF  PRESSURE. 
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VII.  COMPLETE  CONFIGURATIONS 


VII. 1 GENERAL  CONSIDERATIONS 

The  flow  field  about  a wing-body  combination  in  a 
fluid  stream  is  characterized  by  strong  interactions 
between  the  body  and  the  wings.  Although  the  gross  fea- 
tures of  the  flows  around  the  separate  components  are 
largely  the  same  as  described  in  Sections  V.l  and  VI. 1, 
significant  distortions  are  introduced  by  mutual  inter- 
ferences. These  discussions  of  complete  configuration 
behavior  should  be  read  after  those  dealing  with  isolated 
bodies  and  tails.  Separate  discussions  will  be  presented 
for  the  cases  of  low  and  high  angles  of  attack  respectively, 
beginning  with  the  former. 

At  low  angles  of  attack,  the  body  base  usually 
exhibits  a thickened  or  separated  axial  boundary  layer. 

A tail  fixed  to  the  base  will  be  at  least  partially 
immersed  in  this  highly-viscous , low  energy  flow.  As  a 
result,  regions  near  the  tail  root  will  be  subjected  to 
reduced  loading.  It  would  be  expected  then  that,  relative 
to  the  isolated  tail  case,  tail  normal  force  would  be 
reduced,  while  the  spanwise  center  of  pressure  would  move 
outboard.  Since  boundary  layer  thickness  increases  towards 
the  rear  of  the  tail,  this  region  would  be  less  highly 
loaded  than  the  forward  portion  and  the  chordwise  center 
of  pressure  would  move  upstream.  All  of  these  effects 
are  found  in  the  data  to  some  extent. 

The  other  major  class  of  effects  at  low  angles  is 
usually  taken  to  be  due  to  the  influence  of  "potential" 
flow  upwash  around  the  body.  This  often  results  in  an 
effectively  increased  angle  of  attack  on  the  tails,  with 
correspondingly  increased  loading  and  an  inward  shift  of 
the  spanwise  center  of  pressure.  Although  these  effects 
are  probably  present  in  the  low-angle  data,  the  dominant 
influence  seems  to  be  that  of  the  boundary  layer,  as  dis- 
cussed above. 

At  higher  angles,  separated  crossflow  boundary  layer 
effects  become  important.  Not  only  are  there  interactions 
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between  the  lee-side  vortices  from  the  tails  and  the  body: 
the  presence  of  the  tails  "blanks  off"  a portion  of  the 
body  crossflow,  leading  to  further  complicated  interactive 
effects.  The  amplification  of  tail  normal  force  becomes 
reduced  by  downwash  due  to  the  body  vortices  and  complex 
flows  are  induced  near  the  tall -body  junction  by  separating 
and  reattaching  boundary  layers.  Figure  51  shows  an  oil- 
flow  picture  of  the  tail  region  on  the  lee-side  of  the  con- 
figuration at  15°  angle  of  attack.  The  complexity  of  the 
flow  field  is  clearly  shown.  Note  the  downward  direction 
of  the  streamlines  on  the  body,  caused  by  the  lee-side 
vortex  flow.  The  tail  lee-side  vortices  can  also  be  seen 
as  can  the  complex  interactive  region  near  the  tail  root. 
This  picture  shows  clearly  the  difficulties  in  attempting 
a theoretical  treatment  of  such  flows  and  the  need  for  a 
semi-empirical  approach. 


VI X. 2 BODY  CARRYOVER  NORMAL-FORCE  COEFFICIENT 

This  method  permits  calculation  of  the  carryover 
loading  on  a body  due  to  the  presence  of  the  lifting  tails. 
The  method  is  valid  for  angles  of  attack  up  to  15°,  and 
is  restricted  to  body-tail  combinations  at  0°  sideslip. 
Application  consists  of  first  determining  the  carryover 
loading  of  a reference  configuration  and  then  correcting 
this  value  for  effects  of  fin  span,  taper  ratio  and  fin 
position. 

When  a body  has  tail  panels  fixed  on  either  side,  the 
lift  on  these  panels  induces  a certain  amount  of  lift  on 
the  body  region  lying  between  them.  This  lift,  which  is 
usually  termed  "carryover  loading",  can  be  estimated  by 
linear  theory  (ref.  22).  However,  since  the  tail  panels 
of  the  present  configurations  are  located  on  a region  of 
the  body  that  is  strongly  influenced  by  viscous  effects, 
linear  theory  was  found  to  be  inadequate  for  prediction 
of  oarryover  loading. 

The  carryover  lift  usually  acts  in  the  same  direction 
as  that  on  the  tails.  A feature  of  the  present  data  was 
the  existence  of  carryover  loading  acting  in  both  the  same 
and  opposite  direction  to  the  tail  lift  for  the  one  caliber 
base.  In  the  absence  of  a detailed  examination  of  the 
nature  of  the  carryover  loading,  which  is  beyond  the  scope 
of  the  present  investigation,  it  is  difficult  to  describe 
accurately  the  causes  of  the  negative  loading.  The  results 
of  the  flow  visualisation  investigation  did  reveal  the 
presence  of  large  vortices  on  the  one  caliber  base  (which 
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showed  negative  loading)  that  were  not  evident  on  the  two 
or  three  caliber  bases.  This  could  have  led  to  the  effect. 
In  addition,  the  tail  panels  could  have  effectively  blanked 
off  a portion  of  the  crossflow  around  the  body,  thus 
causing  a viscous  load  reduction  which  would  then  appear 
as  a negative  carryover  loading. 

The  carryover  loading  on  the  two  and  three  caliber 
bases  was  essentially  linear  with  angle  of  attack  over  the 
entire  a range  {0°  to  15°) . This  is  consistent  with  the 
solution  from  linear  theory.  However,  because  of  the  flow 
characteristics  over  the  conical  bases,  the  amount  of 
carryover  loading  was  found  to  vary  with  the  location  of 
the  fin  on  the  base,  the  fin  taper  ratio,  and  the  fin  span. 
Therefore,  it  was  decided  to  base  the  method  on  a correla- 
tion of  actual  data  to  describe  a function  of  the  form 


- C (ot,r/s,X,l  p/d) 

NB(T)  NB(T)  LE 


(38) 


Description  of  Method 

It  1b  assumed  that  the  effects  of  the  various  param- 
eters in  equation  (38)  are  separable.  This  allows  the 
equation  to  be  rewritten  as 


'N 


B(T) 


a 

CN 

L NB(T) 


<rQ/s) 


f2(X) 


s ^ *Lh/^ 


(39) 


where  [ CM  is  the  carryover-normal-force-curve  for  a 

l B(T)J  O 

reference  configuration.  The  values  of  the  parameters  for 
the  reference  condition  are 


where 

fck 


r 

•«  o 


N 


B (T)  J 


rQ/s  - 0.583 
X - 0.5 
lLE/d  - 0.00 

is  the  maximum  radius  of  the  body.  The  value  of 
is  given  in  Figure  52  for  each  of  the  three 


bases  tested. 


..  i 
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FIGURE  52.  BASIC  TERM  FOR  C 
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To  determine  values  for  conditions  other  than  the 
reference  conditions  above,  CM  obtained  from  Figure  52 

nB(T5 

should  be  multiplied  by  the  multiplicative  correction 
factors  for  (1)  effect  of  fin  span  from  Figure  53 , (2) 
effect  of  tail  fin  taper  ratio  from  Figure  54,  and  (3) 
effect  of  fin  position  from  Figure  55.  In  determining  the 
ratio  of  body  radius  to  fin  semispan,  r /s,  the  maximum 
body  radius  should  be  used.  J 

Use  of  Method 


To  apply  this  method,  the  following  values  are  required: 
Fins 


Exposed  semispan,  b/2 
Taper  ratio,  X 
Body 


Maximum  body  radius,  rQ 
Base  length,  lb/d 
Fin-Body 


Axial  distance  of  fin  leading  edge  from  start  of 
baBe  section,  lLB/d 

The  procedure  to  be  followed  is: 


1)  Given,  the  base  length,  1^/d,  determine  the  basic 


term  CM 

L NB(T)Jo 
Figure  52. 


from  the  appropriate  curve  of 


2)  Determine  r0/s  ■ ro^ro  + b/2)  then  find  the 
multiplicative  correction  factor  for  fin  span  from 
Figure  53. 


3)  Find  the  correction  factor  for  fin  taper  ratio 
from  Figure  54. 

4)  Given  Vsj/d,  determine  the  correction  factor  for 
the  effect  of  fin  position. 


5)  Determine  CM  from  equation  (39). 

B(T) 
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No  numerical  example  will  be  given  at  this  stage. 
This  will  be  postponed  until  a complete  configuration  is 
considered  in  Section  VII. 7. 


VII. 3 BODY  CARRYOVER  LOAD  CENTER  OF  PRESSURE 


The  following  method  is  used  to  estimate  the  center  of 
pressure  of  the  body  carryover  normal-force  coefficient. 
This  method  is  valid  for  models  of  the  present  type  at 
angles  of  attack  up  to  fifteen  degrees.  Examination  of 
the  data  showed  that,  to  a sufficient  degree  of  accuracy, 
the  centroid  of  the  body  planform  area  between  the  fins 
provided  a suitable  estimation  of  the  center  of  pressure 


As  discussed  in  Section  VII. 2,  the  lift  on  the  tail 
panels  induces  a lift  on  the  body  region  lying  between 
these  panels.  This  interference  lift  is  termed  "carryover 
loading".  The  center  of  pressure  of  this  load  then  lies 
somewhere  between  the  tail  panels,  i.e.  in  the  shaded 
region  of  the  body  in  Figure  56 . 

Experimentally,  the  center  of  pressure  of  CM  is 

NJB  (T) 

found  by  subtracting  the  moments  due  to  the  body-alone  and 
the  tail  panels  from  the  total  moment  and  dividing  by  the 
carryover  normal  force,  as  in  equation  (40)  below. 


Because  of  the  method  required  to  obtain  Xn#  there  can  be 
quite  a lot  of  scatter  in  the  data.  Thus,  it  was  not 
possible  to  maxe  a detailed  examination  of  the  effects  of 
various  parameters,  such  as  fin  taper  ratio,  etc.,  on  XN. 

Using  a lifting-line  potential  flow  model,  Pitts,  at  al. 
(ref.  2)  show  that  5L  varied  with  aepect  ratio,  taper 
ratio,  and  the  ratio  of  body  radius  to  fin  eemiapan.  Due 
to  data  scatter  it  was  not  possible  to  check  the  application 
of  that  model  to  the  present  configurations.  However,  a 
convenient  engineering  method  was  derived  that  does  provide 
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a suitably  accurate  estimate  of  X.,  for  preliminary  design 
purposes.  This  method  is  baaed  onNthe  assumption  that  the 
carryover  loading  on  the  body  is  uniformly  distributed  over 
the  planform  area  of  the  body  region  lying  between  the  tail 
panels.  The  center  of  pressure  of  the  carryover  loading 
is  then  simply  the  centroid  of  that  planform  area. 

Description  of  Method 

For  a conical  base,  the  planform  area  of  the  region  of 
carryover  loading  is  a trapezoid,  as  shown  in  Figure  56. 

The  centroid  of  this  region  is  given  by 


AX 

T 


2r , + r. 


i2  + r, 


(41) 


where  AX  is  the  axial  length  of  the  region,  and  1 is 
measured  from  the  upstream  end  of  the  region. 

Use  of  Method 

To  apply  this  method,  the  body  radii  at  the  fin  root- 
chord  leading  and  trailing  edges,  and  the  axial  extent  of 
the  root  chord  are  required.  These  are  then  substituted 
into  equation  (4.1)  to  determine  lc,  and  henoe  Xm.  A 
numerical  example  will  be  postponed  until  a complete  con- 
figuration is  considered  in  Section  VII. 7. 


VII. 4 TAIL  NORMAL-FORCE  AMPLIFICATION  FACTOR 


This  method  permits  calculation  of  the  amplification 
of  tail  normal  force  when  in  the  presence  of  the  body. 

The  method  is  valid  for  body-tail  combinations  of  the  pres- 
ent type  in  incompressible  flow  and  for  angles  of  attack  up 
to  15°.  From  zero  to  5 the  amplification  is  aonstant. 
Beyond  5 it  increases  linearly. 


The  normal  force  on  a horizontal  fin  fixed  to  a body 
is  usually  greater  than  that  on  an  isolated  fin  at  the  same 
angle  of  attack  due  to  the  upwash  around  the  body.  Hence, 
the  tail  normal-force  amplification,  RT{B) > expressed  sb 

CN, 


^T  (B) 


T(B) 

* 

”N_ 


(42) 
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is  expected  to  be  greater  than  unity.  For  fins  mounted  on 
a body  section  of  constant  diameter,  potential  flow  theory 
(see  ref.  2)  shows  that  the  amplification  is  a function  of 
the  ratio  of  fin  span  to  body  radius  and  varies  between  one 
and  two. 

However,  this  concept  assumes  that  the  flow  over  the 
body  has  not  separated.  When  this  occurs,  the  value  of 

will  depend  upon  how  much  of  the  fin  is  in  separated 

flow.  As  shown  earlier,  the  flow  over  the  bodies  in  the 
present  teat  did  separate  on  the  base  section.  The  result 
was  that  the  experimental  tail  normal-force  amplification 
varied  between  values  greater  than  one  all  the  way  down  to 
zero. 


In  the  latter  case  the  fins  are  in  what  is  essentially 
a completely  separated  flow  and  experience  little  or  no 
force.  This  condition  was  often  found  for  base  B1  below  an 
angle  of  attack  of  5°.  In  general,  is  constant  from 

zero  to  5°  and  then  increases  linearly.  Thus,  the  amplifi- 
cation is  of  the  form 


0 < a < 5° 


a < 15° 


(43) 


When  compared  to  potential  theory , the  present  values 
were  always  less  due  to  a portion  of  the  fins  being  in 
separated  flow  as  described  above.  As  a measure  of  com- 
parison, the  potential-flow  relation  (ref.  2)  for  fins 
mounted  on  a body  of  constant  cross  section  can  be  used  to 
determine  the  radius  of  such  a body  that  would  give  the 
same  amplification  as  in  the  present  test.  Since  this 
relation  gives  a minimum  amplification  of  1.0,  the  equi- 
valent slender  body  radius  can  be  found  only  for  experi- 
mental values  of  Rrn/B\  greater  than  one.  For  these  values, 
the  equivalent  radiusBwas  found  to  be  10-20  percent  of  the 
radius  at  the  root-chord  leading  edge  of  the  fins , and  was 
approximately  half  the  radius  at  the  fin  centroid.  Thus, 
the  tail  normal-force  amplification,  for  body-tail  con- 
figurations used  in  the  present  test,  is  much  less  than 
predicted  by  linear  flow  theory  due  to  flow  separation. 
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The  parameters  affecting  R^,^  include  fin  span/body 

diameter  ratio , aspect  ratio,  taper  ratio,  Reynolds  number 
based  on  maximum  body  diameter,  base  length,  and  the  axial 
location  of  the  fin  on  the  base,  in  addition  to  the  pre- 
viously discussed  dependence  on  the  position  of  separation 
and  the  angle  of  attack.  The  effect  of  axial  fin  location 
of  separation  with  respect  to  the  fin  and  the  effect  of  a 
change  in  body  radius.  Thus,  the  functional  relationship 
for  1^(3)  is  written 

*T(B>  ' ”t(b)  (p*> i-Re-  T’TT'“)  144 

The  method  described  in  the  next  section  was  based  on  the 
data  from  the  present  test. 

Description  of  Method 

Based  on  evaluation  of  the  data,  the  value  of  R-.q.  is 

separated  into  two  terms.  The  first  is  a constant  and  is 

the  value  of  from  zero  to  5°.  This  relationship  is 

expressed  quantitatively  in  equation  (43) . Each  of  the  two 

terms,  Rq  and  Rt , are  determined  by  adding  the  incremental 

effects  of  aspect  ratio,  taper  ratio,  Reynolds  number,  and 

fin  position,  all  for  the  particular  base  under  study,  to 

the  value  at  a given  span  at  given  reference  conditions. 

Thus,  for  instance,  RA  can  be  written  as 

o 


<Ro}d/b 


A*Ro  4 


*XRo  + 


*ReRo  4 


*lRo 


(45) 


The  specified  reference  conditions,  at  which  the  first 
term  was  obtained,  are 


PR  • 1.0 


X - 0.5 


Red 

XLE 


T 


2 .3*10* 

0.0 


The  first  term  is  determined  from  figure  57,  for  the 
particular  values  of  fin  span  and  base  length.  To  this  are 
added  the  effects  of  the  other  parameters,  determined  from 
figures  58  through  61.  In  figure  58,  the  small  vertical 
mark  at  the  end  of  each  curve  of  constant  l^/d  indicates 
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the  limiting  aspect  ratio  for  the  reference  conditions. 
Also,  in  figure  61,  for  the  effect  of  fin  position,  the 
dashed  line  indicates  that,  for  the  given  reference  con- 
ditions, there  is  a maximum  value  of  l-E/d,  beyond  which 
the  fin  is  no  longer  on  the  base. 

The  value  of  Rt  is  determined  in  exactly  the 
manner,  starting  with  the  initial  term  for  the  effect  of 
fin  span  given  in  figure  62.  The  additional  effects  are 
determined  from  figures  63  through  66. 

Use  of  Method 


This  section  demonstrates  the  use  of  the  method  and 
makes  comparisons  with  experimental  data.  The  required 
physical  characteristics  of  the  body-tail  combination  are 
the  following r 

Fins 


Aspect  ratio,  PR 
Taper  ratio,  X 
Span,  b 
Body 

Base  length,  1^ 

Maximum  body  diameter,  d 

Reynolds  number  based  on  maximum  body  diameter,  Red 
Fin- Body 

Position  of  fin  root-chord  leading  edge  from  start 
of  base  section,  1LE 

The  procedure  to  be  followed  isi 

(1)  Given  the  fin  span  and  base  length,  determine  (R0)d/h 

from  the  appropriate  curve  of  figure  57.  ' 

(2)  Determine  the  incremental  effects  of  aspect  ratio, 

taper  ratio,  Reynolds  number,  and  fin  position  from 
figures  58,  59,  60  and  61,  respectively. 

(3)  Use  equation  (42)  to  determine  RQ. 


130 


P/  I SO'O 


NCSC  TH-' 30-78 


FIGURE  62.  EFFECT  OF  FIN  SPAN  ON  SLOPE. 
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EFFECT  OF  ASPECT  RATIO  ON  R, . 
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FIGURE  64.  EFFECT  OF  TAPER  RATIO  ON 
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(4)  For  the  given  values  of  fin  span  and  base  length,  find 

the  value  of  (Ri)  . from  the  appropriate  curve  of 
figure  62.  a/D 

(5)  As  in  step  2,  find  the  incremental  effects  of  aspect 

ratio,  taper  ratio,  Reynolds  number,  and  fin  position 
from  figures  63,  64,  65,  and  66,  respectively. 

(6)  The  slope,  Rj,  is  determined  from  the  equation 


R.  - (R.»d/b  + 4*R1  + 4Ari  + 4ReR>  + 41R.  (46> 

(7)  The  value  of  R^^j  is  given  by  aquation  (43),  repeated 
here . 


^(B)  “ { 


0 1 a <L  S' 


R + Rj  (a)  , 5°  < a £ 15° 


Numerical  Example.  This  example  compares  the  method  to 
data  from  the  present  test.  The  configuration  chosen  is 
N2C2B2T9m  at  a free-stream  Reynolds  number  of  2.3x10',  based 
on  the  maximum  body  diameter.  The  values  of  the  required 
parameters  are: 


Fins 

Aspect  ratio i 
Taper  ratio: 
Span: 

Body 


PR  - 2.0 
A » 0 . 0 
b - 5.0  in. 


Bate  length:  lfa  ■ 14.0  in. 

Maximum  body  diameter:  d « 7.0  in. 

Reynolds  number  based  on  maximum  body  diameter: 
Red  - 2.3x10* 


NCSC  TM-238-78 


4 


Pin  Body 

Position  of  fin  root-chord  leading  edge  from  start  of 
base  section:  1Le  - 0.0  inch.  A sketch  of  this  configura-  [ 

tion  is  given  in  figure  67. 

(1)  Prom  figure  57,  (R^a/iy  “ 0*927. 

(2)  From  figures  58-61  the  following  values  are  obtained: 


4*Ro  - 

0.37 

from  figure  58 

■ ■ 

0.18 

from  figure  59 

VRo  ' 

0.00 

from  figure  60 

41R0  ■ 

0.00 

from  figure  61 

(3) 

Using  equation 

. (43)  , 

Rq  - 1.1170 

(4) 

Prom  figure  62 

' # (Rj ) 

4/b  " °-0135 

(5) 

Using  figures 

63-66, 

the  following  values  are  determined 

A*Rl  * “ 

0.016 

from  figure  63 

AXRl  - 

0.017 

from  figure  64 

AReRi  - 

0.00 

from  figure  65 

- 

AlRl  " 

0.00 

from  figure  66 

(6)  Using  equation  (46),  R,  has  the  value  0.0136  deg"1. 
Thus,  the  value  of  is 
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FIGURE  68.  TAIL  NORMAL-FORCE  AMPLIFICATION  COMPARISON 
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VII. 5 TAIL-ON-BODY  CHORDWISE  CENTER  OF  PRESSURE 

This  method  permits  estimation  of  the  chordwise  loca- 
tion of  center  of  pressure  on  tail  panels  when  mounted 
horizontally  on  bodies  with  conical  bases.  The  method  is 
valid  for  low-aspect-ratio  panels  with  rounded  leading 
edges  and  blunt  trailing  edges  at  body  angles  of  attack 
up  to  15°  in  incompressible  flow.  It  is  developed  as  a 
simple  correction  to  the  chordwise  center  of  pressure  on 
an  isolated  tail  panel  section  VI. 3.  It  is  recommended  that 
the  method  be  used  only  for  tails  on  conical  bases  of 
slenderness  ratios  greater  than  1 caliber.  A separate 
discussion  is  presented  concerning  the  utility  of  tail 
panels  mounted  on  a cne  caliber  base. 

As  discussed  in  Section  VI. 3,  a correlation  of  actual 
data  for  the  center  of  pressure  on  isolated  tails  is  a more 
satisfactory  method  than  the  various  linear  theories.  This 
is  found  to  be  the  case  also  for  tails  in  the  presence  of 
the  body.  Results  from  linear  theory,  as  described  in 
reference  2,  indicate  very  little  difference  in  the  loca- 
tion of  the  chordwise  center  of  pressure  for  the  isolated 
tail  and  the  tail  in  the  presence  of  the  body.  However, 
the  data  from  the  present  test  shows  that  is 

slightly  forward  of  Jfij,  the  difference  decreasing  with 
increasing  angle  of  attack.  In  previous  sections  the 
possibility  and  consequences  of  the  oxistenoe  of  a thick 
boundary  layer  or  a region  of  separated  flow  over  the  aft  ; 

portion  of  the  base  near  zero  angle  of  attack  has  been 
discussed.  If  this  type  of  flow  was  present,  then  the  aft 
portion  of  the  tail  panels  would  be  embedded  in  this  sepa-  ' 

rated  flow.  Then,  the  rear  portion  of  the  tail  panels  * 

would  generate  reduced  lift  and  the  center  of  pressure 
would  move  forward,  consistent  with  the  present  data. 

Examination  of  the  data  revealed  a consistent  forward  j 

shift  of  X_ from  the  isolated  tail  data.  This  shift 
was  a function' of  the  angle  of  attack  but  was  almost  inde-  . 

pendent  of  model  configuration.  Hence,  an  accurate  estimate 
of  5? 
tion 

this  is  limited  to  tails  mounted  on  the  two  and  three 
caliber  bases  only. 

Chordwise  Center  of  Pressure  of  Tails  Mounted  on  the 
fine  da  liber base""  "" 


could  be  determined  by  applying  a simple  correo- 
to  the  estimate  of  Xm  for  the  same  tail.  However, 


The  data  for  tail  panels  mounted  on  the  one  caliber 
base  showed  a tremendous  amount  of  scatter,  as  can  be  seen 
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in  figure  69.  This  prevented  any  correlation  of  the  one 
caliber  base  data.  Referring  to  the  oil  flow  photograph 
in  Section  V.l,  it  is  possible  to  understand  the  cause  of 
this  catter.  The  large  vortices  on  the  base  greatly  affec- 
ted the  loads,  both  on  the  base  and  on  the  tail  panels. 

As  was  shown  in  Section  VI l. 2,  the  result  could  be  the 
almost  total  elimination  of  the  tail  normal  force,  indica- 
ting that  the  tail  panels  were  completely  embedded  in 
separated  flow.  This  does  not  provide  a very  desirable 
situation  for  vehicle  control.  Therefore,  it  is  recommended 
that  the  one  caliber  base  not  be  used  unless  specifically 
required.  If  its  use  is  mandated  then  data  from  this  test 
(not  given  here)  should  be  used  to  animate  XT^gj. 

Description  of  Method 

The  method  involves  applying  a simple  correction  to  the 
estimate  for  the  tail-alone  chordwise  center  of  pressure, 

$p,  for  the  same  tail  panel.  Thus, 

XT(B)  - XT  - AS f (48) 

The  variation  of  a5c  with  angle  of  attack  is  given  in 
Figure  70.  The  proper  value  of  3?T  should  be  determined 
by  the  method  described  in  Section  VI. 3. 

Use  of  Method 


This  section  demonstrates  the  use  of  the  method  and 
makes  comparisons  with  experimental  data.  The  required 
physical  parameters  are  the  following » 

Taper  Ratio 

Aspect  Ratio  (baaed  on  complete  symmetric  planform) 

(1)  Determine  the  variation  of  XT  with  angle  of  attack 
from  Section  VI. 3. 

(2)  Apply  the  correction  given  in  Figure  70  to  determine 
XT(B),  using  equation  (48). 

Numerical  Example.  This  example  compares  the  method 
to  data  from  the  present  test.  The  configuration  chosen 
is  N2C2B2Tl2u  at  a free-stream  Reynolds  number  of  16.1x10*. 
The  values  of  the  required  parameters  are 


Taper  Ratio  t X - 0 . 5 
Aspect  Ratio » PR  ■ 2.0 
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Using  Figure  47(c)  of  Section  VI. 3 and  Figure  70  of  the 
present  section  the  following  values  are  obtained. 


o£ 

VCR 

a*/cR 

0 

.453 

.08 

.373 

5 

.435 

.06 

.375 

10 

.430 

.04 

.390 

15 

.435 

.02 

.415 

These  values  are  compared  with  experiment  in  Figure  71. 
The  agreement  is  seen  to  be  quite  satisfactory. 


VII. 6 TAIL-ON -BODY  SPANWISE  CENTER  OF  PRESSURE 

This  method  permits  estimation  of  the  spanwise  loca- 
tion of  center  of  pressure  on  tall  panels  mounted  horizon- 
tally on  bodies  with  conical  bases.  The  method  is  based 
om  strip  theory  and  assumes  a constant  value  of  the  section 
normal-force-curve  slope  along  the  span.  It  is  found  to 
be  valid  throughout  the  range  of  experimental  values  tasted. 

In  contrast  to  the  chordwise  location  of  center  of 
pressure  (Section  VII. 5),  the  spanwise  location  is  expeoted 
to  be  different  for  a tail  panel  in  the  presence  of  the 
body  than  for  an  isolated  tail.  The  inner  portions  of  the 
panels  are  in  the  accelerated  flow  around  the  body  whereas 
the  outer  portions  are  in  a flow  that  approaches  free-stream 
conditions.  Hence,  the  flow  speed  normal  to  the  panel  is 
not  uniform  along  the  span,  being  greater  near  the  root 
chord.  This  results  in  an  effectively  greater  angle  of 
attack,  and  therefore,  a correspondingly  greater  load,  on 
the  inner  portion  of  the  panel.  The  final  result  is  that 
the  spanwise  location  of  center  of  pressure  is  inboard  on 
a wing-body  combination  compared  to  an  isolated  wing.  As 
the  size  of  the  body  increases  with  respect  to  the  tail 
semispan  (increasing  r/s) , than,  eventually,  the  entire 
panel  becomes  subjected  to  an  essentially  constant  flow 
field  and  the  spanwise  center  of  pressure  moves  outboard 
again. 

'’his  concept  has  been  expressed  quantitatively  in 
reference  23.  Using  strip  theory,  ignoring  thickness 
effects,  and  assuming  a constant  section  normal-force- 


144 


*T(B) 


n.  *«Hiiibiujimi3v 


• . 1 3-^1  (vw  « navyoys^  iftfrMwv*  Wu’fW't  - 


NCSC  TM-238-78 


0.08 


riOUWt  71,  COMPARISON  BETWEEN  PREDICTED  AND  EXPERIMENTAL 
CHORDWISE  CENTER  OP  PRESSURE  OP  TAIL 


NCSC  TM-238-78 


curve 

whose 

total 


where 

Cy  is  the  chord  length , and  r is  the  body  radius.  The 
expression  for  au  is 

au  - (1  + r2/y* ) a (50) 

for  an  unrolled  body,  where  a is  the  angle  of  attack  of 
the  body.  Similarly,  the  rolling  moment  can  be  expressed 

as 

s 

RM  - q cn  f ou  Cy(y  - r)dy  (51) 

r 


slope,  Cn  , equal  to  that  of  a two-dimensional  wing 

a 

sweep  is  the  same  as  that  of  the  elemental  strip,  the 
normal  force  on  the  lifting  surface  was  found  to  be 

8 


NP 


q c. 


n 


oi  * 


<XU  Cy  dy 


(49) 


q is  the  dynamic  pressure,  ctu  is  the  upwash  angle, 


Carrying  through  the  integrations  and  dividing  the  rolling 
moment  by  the  normal  force,  the  following  expression  is 
obtained  for  the  spanwise  canter  of  pressure,  measured 
from  the  root  chord. 
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where  A is  the  fin  taper  ratio  and  b/2  is  the  semispan 
of  the  exposed  fin.  This  relation  is  shown  graphically  in 
Figure  72,  taken  from  reference  23,  for  three  values  of  A. 

Description  of  Method 

When  equation  (52)  was  compared  with  data  from  the 
present  test  the  agreement  was  found  to  be  quite  satisfac- 
tory. In  reference  23  agreement  with  data  in  transonic 
flow  was  also  seen  to  be  very  satisfactory.  Thus,  this 
theory  seems  to  apply  over  a wide  range  of  flow  conditions. 


In  the  case  of  tails  mounted  on  conical  bases,  the 
appropriate  value  of  body  radius  must  he  determined.  A 
brief  study  of  various  radii  (such  as  at  the  mid  point  of 
the  root  chord)  was  conducted.  However,  Figure  72  shows 
that  is  relatively  insensitive  to  the  body  radius. 

This  means' that  when  applied  to  bodies  of  the  present  type, 
the  most  convenient  radius  can  be  used.  The  radius  at  the 
leading  edge  of  the  root  chord  is  very  convenient  to  deter- 
mine and  yields  quite  good  results.  Thus,  this  radius  is 
recommended  in  the  use  of  this  method. 


Use  of  Method 


This  section  demonstrates  the  use  of  the  method  and 
makes  comparison  with  experimental  data.  The  required 
physical  parameters  are  the  following: 

Body  radius  at  root-chord  leading  edge  of  fin,  rLE 

Semi span  of  fin,  s 

Fin  taper  ratio,  A 


(1)  Determine  ratio  of  body  radius  at  root-chord  leading 
edge  of  fin  to  fin  semiapan,  r/s. 

(2)  Find  as  a fraction  of  exposed  semispan  from 

the  appropriate  curve  of  Figure  72. 


Numerical  Example.  This  example  compares  the  method 
to  data  from  the  present  test.  The  configuration  chosen 
is  N2C2B2T8A  at  a free-stream  Reynolds  number  of  16. lx  10 8 . 
The  values  of  the  required  parameters  are: 


t 


Body  radius  at  root-chord  leading  edge  of  fin: 

'Ll  - 3-S  ln- 
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Semispan  of  fin:  s = 6.0  in. 

Taper  ratio:  X « 0.0 

(1)  The  body  radius  to  fin  semispan  ratio  is 

rLE/s  " 0.583 

(2)  Froj|  the  lower  curve  of  Figure  72  (X  ■ 0)  , the  value 
of  YT(B)/b/2  is  0.35. 

This  value  ie  compared  with  the  experimental  values  in 
Figure  73.  Agreement  is  quite  good.  Additional  compari- 
sons are  given  in  Figure  74  for  the  configuration  N2C2B2T12 
with  fins  mounted  in  the  forward  (Figure  74a),  mid  (74b), 
and  aft  (74c)  positions.  The  agreement  is  seen  to  be  very 
good  in  all  oases. 


VII. 7 USE  OF  METHODS  TO  PREDICT  COMPLETE  CONFIGURATION 
CHARACTERISTICS 

In  this  section  the  previously  developed  methods  are 
aombined  to  determine  the  total  body-tail  normal-force 
coefficient  and  pitching-moment  coefficient  for  several 
complete  configurations.  The  individual  methods  are  worked 
out  in  detail  for  one  configuration  which  then  serves  as 
an  overall  example  for  analysis  of  others. 

Detailed  Numerical  Example 

The  configuration  chosen  for  detailed  analysis  is 
N2C2B2T11{4  at  a free-stream  Reynolds  number  of  16.38*10  . 
The  components  of  this  configuration  are,  in  a sense,  a 
geometric  average  of  all  the  components  tested. 

Although  only  the  total  normal  force  coefficient  and 
pitching  moment  coefficient  will  be  compared  with  data  in 
this  section,  all  forces  and  moments  for  which  a method 
exists,  will  be  oomputed  in  detail.  This  then  serves  as 
a complete  example  of  applying  all  the  methods  developed 
in  the  previous  seations. 

Values  of  Required  Parameters 

Fin 


« - 1.0 
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Pin-Body 


1LE'/d 


0.00 


Application  of  Methods 


'N 


T 


(Section  VI. 2) 

Prom  figure  40 
Prom  figure  41 


C„  - 1.570  rad 
Na 

a,  - 0.657  rad 


-» 


-2 


The  tail-alone  normal-force  coefficient  is  now  deter- 
mined from  equation  (33) , given  below. 


X 


- C, 


'N. 


• a + a. 


T "a 

The  results  are  tabulated  in  Table  II. 

XT  (Section  VI. 3) 

The  tail-alone  chordwise  center  of  pressure  is  deter- 
mined from  figure  47(b)  and  the  results  are  tabulated  in 
Table  II. 


Nb 

(1) 

(2) 


(3) 

(4) 

IS) 

(6) 


(Section  V. 2) 

From  figure  17,  xBa/ljj  ■ 0.455. 

The  equation  for  the  diameter  of  a cone  two  calibers 
in  length  is 

dx  - d (1  - |) 

where  d is  the  diameter  of  the  base  of  the  cone. 
Thus  db  - 3.815  in.  and  dfa/d  - 0,545. 

From  equation  (13),  CM  ■ 0.594  per  radian. 


Use  C, 


N. 


Na 

■ Cjj  a to  find  potential  force  as  a function 


pot  "a 
of  angle  of  attack. 

Determine  x from  equation  (14)  for  the  angles  of 
attack  Xieteacin  Table  II. 

Calculate  the  planform  axes  from  each  x in  step 
(5)  to  the  end  of  the  body. 
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(7)  Find  the  appropriate  value  of  C.  for  each  local 

Qo 

section  at  each  crossflow  Reynolds  number  from  figure 
18.  In  this  example  Cd  - 0.290  throughout. 

(8)  Use  equation  (15)  to  find  CM  . 

Nr.l 

(9)  Add  result  from  step  (4)  to  result  from  step  (8)  to 
find  CM  . 

nB 

Since  this  configuration  is  the  one  used  as  the  detailed 
numerical  example  in  Section  V. 2,  only  the  final  values  of 
CM  are  listed  in  Table  II. 
nb 

CM  (Section  V. 3) 


(1)  From  figure  25,  (Ka  - KJ  - 0.898 

(2)  From  figure  26,  Km  - 0.940 

(3)  Using  figure  17  of  Section  V.2,  x /l.  - 0.455.  Since 
lb/d  - 2 for  base  B2,  xga/d  - 0.9Io,  measured  from  the 

start  of  the  base  section.  Measured  from  the  nose, 
xmm /d  ■ 5.91  -*■  x_.  - 3.4475  ft,  This  is  the  upper 

limit  of  integration  in  equation  (26) . 

(4)  Substituting  the  necessary  quantities  into  equation 
(26)  and  carrying  out  the  integration 

- 8.043  rad"1 


(5) 

From  figure 

27,  n - 0.635 

(6) 

Find  xgc/d 

at  7°,  10°,  13°,  from  equation  (28) 

of, 

xsc/d 

7 

2.667 

10 

1.833 

13 

1.556 
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(7) 

The  crosaf low-drag  coefficient,  from  figure  18  of 

Section  V.2  is  0.290  for  all  sections 

(8) 

The  results  of  the  integration  of  equation  (27)  are 

4* 

7 

-0.4367 

10 

-0.0091 

13 

+ 0.0963 

The  value*  in  Table  II  were  determined  by  linearly  interpola- 
ting out  to  13°  and  assuming  a constant  value  of  CM 
beyond  13®.  Mnl 

CA 

ab 

(Section  V.4) 

(1) 

(2) 

From  figure  31. 

Using  figure  32. 

Cn  - 0.0258 
°b 

Cg  * .00278 

(3) 

From  equation  (31).  C_  - 0.0685 

(4) 

Combining  the  results  of  steps  (1)  and  (3) , 

i 

(CA  ) - C 

aB  o °b 

+ C - .0943 
°fi 

1 

i 

(5) 

From  figure  33. 

a1  ■ . 109  rad"1 . Thus, 

i 

n 

C.  - 0.0943 
*B 

+ .109  u 

j 

The 

results  are  listed  in  Table  II. 

i 

CM  (Section  VII. 2) 

NB(T) 

4 . 

H 

• 

(1) 

From  figure  52, 
1.4209  rad"1 

C„  1 is  linear  with  a slope  of 

nB(T)J  o 

t 1 

\ 

(2) 

r0/s  - 3. 5/(3. 5 

+ 2.5)  - 0.5B3 

*'  t 

1 

From  figure  53, 

f j (re/s)  ■ 1.00 

J 
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(3)  The  correction  factor  for  fin  taper  ratio  from  figure 
54  ia  f2  (A)  - 1.00 

(4)  From  figure  55,  f,(lLE/d)  ■ 1.00 

(5)  Uaing  equation  (39) , 

C.,  ■ 1.4209  a 

nb<t> 


The  values  of  the  carryover  loading  are  tabulated  in  Table 
II. 


^T  (B) 

(Section  VII. 4) 

(1) 

From 

figure  57,  (R0)tf//b  ■ 0.927 

(2) 

From 

figure  58,  - 0.000 

(3) 

From 

figure  59,  A^R0  ■ 0.000 

(4) 

From 

figure  60,  AR#Ro  ■ 0.000 

(5) 

From 

figure  61,  A^RQ  ■ 0.000 

Thus 

R « 0.927 
o 

(6) 

From 

figure  62,  < R i ) dl/b  " 

(7) 

From 

figure  63,  R}  * 0.000 

(8) 

From 

figure  64,  A^Rj  ■ 0.000 

(9) 

From 

figure  65,  Aps|R,  ■ 0.000 

(10) 

From 

figure  66,  A^R,  - C.000 

Thus 

Rj  - 0.0135  deg“l 

Combining 

these  results, 

[ 0.927 

a £ 5' 

*T(B)  “ \ 

0.927  + 0.0135  a 

a > 5 

These  result*  are  listed  in  Table  II. 
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(Section  VII. 3) 

The  body  radius  at  the  root  chord  leading  edge  is  3.5 
in.  To  determine  the  body  radius  at  the  root  chord  trailing 
edge,  first  determine  the  axial  length  of  the  root  chord. 
Referring  to  the  sketch  below 


where  1 is  the  half  angle  of  the  vortex  of  the  cone 
planform. 

For  a two  caliber  base,  y - 14.04°.  Thus 
xa  - 6.665  cos  (14.04) 

" 6.465  in. 

The  equation  for  the  base  radial  distribution  is 


r<«>  - 1 U - 


starting  at  the  beginning  of  the  base  section.  Therefore, 
the  body  radius  at  the  root-chord  trailing  edge  is  1.864  in. 
The  equation  for  the  centroid  of  a trapezoid  is 


-<P 


2ra  + r j 

ra  + ri 


where  r2  - 3.5  in. 

ra  - 1.884  in. 
dx  « 6.465  in. 

Substituting  these  values  into  the  last  equation,  the  axial 
distance  of  the  centroid  from  x^g  is  1Q  ■ 2.909  in. 


j 
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W 


The  axial  distance  between  the  moment  reference  center 
™ 30.5  in.)  and  the  root  chord  leading  edge  (start  of 
base  section)  is  4.5  in.  Hence,  XM/d,  the  distance 
between  the  moment  reference  and  tfle  centroid  is 

5^/d  - — — -+  y2 ,9Q9  - 1.0584 
(Section  VII. 5) 


mu  -Th®re.are  two  •t®P8  to  the  determination  of  Xm/_. . 
The  first  is  the  application  of  the  correction  giveX^f^ 
figure  70  to  Xp.  The  second  is  the  conversion  of  this 
result  to  axial  distances  from  the  moment  reference  center. 
JJ*  Procedure  for  the  first  part  is  outlined  in  Seotion 
VII. 5.  To  determine  the  axial  distance  from  the  moment 
reference  center,  refer  to  the  following  sketch. 


The  equation  for  X,p,_./d,  measured  from  the  moment 
reference  center,  theR' is 


5T(b)  (Section  VII. 6) 

(1)  r/s  - 0.583 

7 

(2)  Prom  figure  70,  ■ 0.420,  independent  of  angle  of 

attaak.  * 
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The  value  of  > is  required  to  determine  the  axial 

diatance  from  tneBmoment  reference  to  and  is  given 

in  Table  II.  ' ' 


Total  Body-Tall  Normal-Force  Coefficient 

The  total  body-tail  normal-force  coefficient,  CM  , 

WBT 

is  given  by  the  following  relation  for  undeflected  fins. 


(53) 


In  making  comparison*  to  the  present  test,  due  to  flow  and 
model  symmetry  it  can  be  assumed  that 


**T  (B)  1 " ^<13)2  " ’SMB)  (54 

Hence 

ST 

CM  - CM  + 2C„  R,p  *=-  + c.  (55 

N3T  NB  NT  P(B)  ®R  ^B(T) 

Combining  the  relevant  values  in  Table  II,  the  variation 
of  CM  with  angle  of  attack  is  obtained.  This  variation 
WBT 

is  tabulated  in  Table  IT  and  is  compared  with  experimental 
results  in  figure  75.  Agreement  between  method  and  experi- 
ment is  very  good  except  at  the  medium  values  of  angle  of 
attaak  where  the  method  over -predicts  CM  . A detailed 

nBT 

examination  of  the  individual  loads  contributing  to  CM 

WBT 

revealed  that  most  of  the  error  was  from  CM  . It  is 

WB 

possible  that  the  body-alone  normal-foroe  coefficient  may 
require  a$  angle  of  attack  term  of  higher  order  than  the 
present  a . Additional  comparisons  are  given  in  figures 
76,  77  and  78.  Thsss  configurations  wsrs  ohossn  to  show 
ths  offsets  of  fin  position  (figurss  75  and  76) , base 
(figures  75  and  77)  , and  fin  span  (figurss  77  sr.d  78) . 

The  agreement  in  figures  76  end  77  is  excellent  but  the 
effect  of  fin  span,  shown  in  figure  78,  is  somewhat 
over-predicted. 
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a{«) 


FIGURE  75.  COMPARISON  BETWEEN  PREDICTED  AND  EXPERIMENTAL 
TOTAL  BODY-TAIL  NORMAL  FORCE  COEFFICIENT. 
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FIGURE  77.  COMPARISON  BETWEEN  PREDICTED  AND  EXPERIMENTAL 
i TOTAL  BODY-TAIL  NORMAL  FORCE  COEFFICIENT. 

i 
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Total  Body-Tall  Pitching-Moment  Coefficient 

The  total  body-tail  pitching-moment  coefficient,  C 
la  given  by  the  following  relation  for  undeflected  fins. 


c ■ c 

«BT  11 


B 


2 


CNTRT(B)| 


XT(B)i  ST  h 

" V*  + C«B(T)  * 


T 


mbt' 


(56) 


Making  the  assumption  that 


'i 

f 

& 


VT  (B)  1 


XT(B)2  “ XT(B) 


(57) 


(58) 


and  combining  equations  (54)  , (56) , and  (57)  yields 

C«BT  ’ % + *V*<«  ^ i + S«,  ^ 

These  values  are  also  tabulated  in  Table  II. 

The  comparison  with  experiment  is  shown  in  figure  79. 
Agreement  is  excellent  throughout  the  angle  of  attack 
range.  Additional  comparisons  are  given  in  Figures  BO, 

81,  and  82  for  the  same  configurations  previously  shown 
in  76-78,  respectively.  The  effect  of  fin  position 

is  only  slightly  over-predicted  in  figure  80  and  the  effect 
of  a different  base  is  handled  well  in  figure  81.  However, 
the  effect  of  fin  span  on  Cu  , figure  82,  appears  to  be 

MBT 

substantially  under-predicted.  By  comparing  figures  81 
and  82,  it  is  seen  that  the  effect  of  larger  fins  is  to 
greatly  reduce  the  value  of  CM  . Thus,  it  appears  possi- 

oT 

ble  that  the  cause  of  the  discrepancy  in  figure  82  is  the 
over-prediction  of  Cu  , the  pitching  moment  due  to  the 

(B) 

tail  in  the  presence  of  the  body.  Referring  to  equation 
(58) , Cu  is  the  second  term  on  the  right  hand  side 

”T  (B) 

and  all  the  components  that  make-up  cu  are  well 

x (B) 

predicted.  Thus,  the  cause  of  the  discrepancy  must  lie  in 
either  Cu  or  Cw  . Examining  the  last  term  first, 

“b  _mb  (T) 

both  CM  and  X.,  are  reasonably  well  predicted. 
nB  (T)  " 
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Therefore,  by  process  of  elimination,  the  discrepancy  is 
due  to  the  prediction  of  C„  . Recalling  the  method  for 

estimating  Cu  , comparisons  with  data  showed  that  the 
mb 

experimental  values  of  the  nonlinear  term  were  slightly 
more  negative  than  the  values  found  by  applying  the  method. 
However,  in  figure  82  just  the  opposite  is  the  case.  The 
methods  for  Cu  and  CM  use  an  affective  separation  point 
mb  nb 

concept.  If  this  affective  separation  point  were  to  move 
downstream  when  the  fins  were  attached  then  there  would  be 
additional  loading  over  the  aft  portion  of  the  body  that 
has  not  bean  accounted  for  in  the  present  methods.  This 
would  increase  the  total  pitahing  moment  over  the  prediction, 
consistent  with  the  data. 
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VIII.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  work  described  in  this  report  demonstrates  the 
utility  and  power  of  the  systematical  experimental  approach 
to  semi-empirical  prediction  method  construction.  Through 
a combination  of  experiment  and  theory,  methods  are  now 
available  for  estimating,  with  good  accuracy,  the  hydro- 
dynamic characteristics  of  certain  types  of  submersible 
vehicles.  At  present  the  methods  are  applicable  only  to 
vehicles  having  undeflected,  cruciform  tails,  at  pitch 
angles  of  attack  to  15  degrees.  The  techniques  can  be 
easily  extended  to  deflected  tails  through  generation  and 
correlation  of  appropriate  data.  The  models  as  fabricated 
will  permit  systematic  tail  deflections.  In  this  way,  the 
return  on  initial  investment  can  be  considerably  increased. 

With  developments  in  vehicle  performance  requirements, 
it  may  become  necessary  to  maneuver  at  combined  pitch  and 
roll,  or  yaw,  or  to  deploy  forward  control  surfaces  to 
enhanae  maneuverability  or  stationkeeping.  The  models, 
at  present,  are  not  flexible  enough  to  permit  testing  of 
such  concepts.  However,  they  can  be  made  to  do  so  in  a 
relatively  straightforward  manner. 

The  most  obvious  required  extension  to  the  developed 
methods  is  the  addition  of  propulsion  effects.  To  accom- 
plish this  would  require  considerable  modification.  It  is 
thought  that  in  view  of  the  configurational  developments 
which  may  be  desirable,  the  addition  of  power  would  best 
be  postponed  until  some  of  the  other  concepts  mentioned 
above  have  been  investigated. 

Recommendations 


It  is  recommended  that  the  following  developments  of 
this  initial  work  be  considered i 

i)  Modify  the  models  such  that  forward  controls  and  com- 
bined pitch-roll  attitudes  can  be  attained.  Mount 
the  forward  controls  on  individual  balances.  Automats 
attitude  and  deflection  changes  to  the  greatest  extent 
possible.  In  this  way,  test  time  will  be  reduced  and 
data  quantity  increased. 
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ii)  Modify  the  models  such  that  power  effects  are  included. 
This  may  be  dona  following  i) , or  immediately.  A 
systematic  investigation  of  power  parameters  will  be 
necessary  for  a wide  range  of  geometric  configurations 
and  it  may  be  thought  advisable  to  postpone  this  step 
until  i)  has  been  completed. 
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